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[Scope of the Claim] 
[Claim 1] 

A tomogram image processing method, comprising the steps of: 

a) carrying out an interpolation operation on a plurality of tomogram images 
of a test boy sliced in an area at predetermined intervals to obtain stereoscopic image 
data elements; 

b) dividing said area into small stereoscopic segments; 

c) counting the number of said stereoscopic image data elements for each of 
said small stereoscopic segments; and 

d) binarizing said small stereoscopic segments on the basis of the number of 
said stereoscopic image data elements. 

[Claim 2] 

A tomogram image processing apparatus, comprising: 

interpolation carrying out means (3) for carrying out an interpolation 
operation on a plurality of tomogram images of a test boy sliced in an area at 
predetermined intervals to obtain stereoscopic image data elements; 

area dividing means (7) for dividing said area into small stereoscopic 
segments; 

stereoscopic image counting means (8) for counting the number of said 
stereoscopic image data elements for each of said small stereoscopic segments; and 

binarizing means (9, 10) for binarizing said small stereoscopic segments on 
the basis of the number of said stereoscopic image data elements. 
[Claim 3] 

An organism magnetic field measuring method, comprising the steps of: 

e) carrying out a plurality of Biot-Savart law calculations on a plurality of 
magnetic field information elements respectively for a plurality of vertexes of said 
small stereoscopic segments obtained in the tomogram image processing method as set 
forth in claim 1; 

f) adding up the results of the Biot-Savart law calculations to obtain a 
calculated magnetic field value; 

g) calculating a difference between said calculated magnetic field value and a 
measured magnetic field value; 

h) correcting data on said magnetic field sources included in the results of the 



1 



Biot-Savart law calculations in accordance with said difference; 

i) judging whether or not said difference is vanishingly small; 

j) repeating said step (e) through step (h) when it is judged in said step (i) that 
said difference is not vanishingly small; and 

k) outputting, as an organism magnetic field-measured outcome, said data on 
said magnetic field sources included in the results of the Biot-Savart law calculations 
carried out in said step (h) when it is judged in said step (i) that said difference is 
vanishingly small. 
[Claim 4] 

An organism magnetic field measuring method as set forth in claim 4, in 

which 

said step (h) is replaced with the steps of: 

hi) calculating a difference between a coordinate of each of said vertexes of 
said small stereoscopic segments and a coordinate of each of corrected vertexes of said 
small stereoscopic segments; 

hi) calculating said coordinate of each of said corrected vertexes of said small 
stereoscopic segments on the basis of a code of said difference calculated in said step 
(hi); and 

h3) redistributing an estimate value of said magnetic field source in 
accordance with said difference calculated in said step (hi) for each of said corrected 
vertexes of said small stereoscopic segments calculated in said step (h2). 
[Claim 5] 

An organism magnetic field measuring apparatus, comprising: 

calculating means (111, 112, 11m; 101, 102, 10m) for carrying out a 
plurality of Biot-Savart law calculations on a plurality of magnetic field information 
elements respectively for a plurality of vertexes of said small stereoscopic segments 
obtained in the tomogram image processing method as set forth in claim 1; 

adding means (20; 200) for adding up the results of the Biot-Savart law 
calculations outputted from said calculating means (111, 112, 11m; 101, 102, 
10m) to obtain a calculated magnetic field value; 

error calculating means (30; 300) for calculating a difference between said 
calculated magnetic field value outputted from said adding means (20; 200) and a 
measured magnetic field value; 

correcting means (111a, 112a, lima) for correcting data on said magnetic 
field sources included the result of the Biot-Savart law calculations outputted from said 
calculating means (111, 112, 11m; 101, 102, 10m) in accordance with said 
difference; and 

corrected result collecting means (40; 500) for collecting the results of the 
Biot-Savart law calculations corrected by said correcting means (111a, 112a, lima) 
and outputting them therethrough as an organism magnetic field-measured outcome. 
[Claim 6] 

An organism magnetic field measuring apparatus as set forth in claim 5, 
further comprising: 
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as correcting means, difference calculating means (41i) for calculating a 
difference between a coordinate of each of said vertexes of said small stereoscopic 
segments and a coordinate of each of corrected vertexes of said small stereoscopic 
segments; 

corrected vertex extracting means (42i) for calculating said coordinate of each 
of said corrected vertexes of said small stereoscopic segments on the basis of a code of 
said difference calculated by said difference calculating means (41i); and 

redistributing means (43i, 44i, 45i) for redistributing an estimate value of said 
magnetic field source in accordance with said difference calculated by said difference 
calculating means (41i) for each of said corrected vertexes of said small stereoscopic 
segments calculated by said corrected vertex extracting means (42i). 

7. An organism magnetic field-measured outcome display method of 
displaying an organism magnetic field-measured outcome along with an organism 
structure reference image, comprising the steps of: 

1) calculating absolute values of said organism magnetic field-measured 
outcomes for a plurality of vivo points; 

m) setting up a size of a reference graphic for each of said vivo points on the 
basis of said absolute values calculated in said step (1); and 

n) displaying said reference graphic whose sizes is set up in said step (m) 
along with said organism structure reference image for each of said vivo points. 

8. An organism magnetic field-measured outcome display method of 
displaying an organism magnetic field-measured outcome as set forth in claim 7, in 
which said reference graphic is in the form of a circular shape. 

9. An organism magnetic field-measured outcome display apparatus for 
displaying organism magnetic field-measured outcomes along with an organism 
structure reference image, comprising: 

reference graphic size calculating means (700) for calculating absolute values 
of said organism magnetic field-measured outcomes for a plurality of vivo points, and 
setting up a size of a reference graphic for each of said vivo points on the basis of said 
absolute values thus calculated; 

organism structure displaying means (800) for displaying said organism 
structure reference image; and 

reference graphic displaying means (900) for displaying said reference 
graphic whose size is set up by said reference graphic size calculating means (700) for 
each of said vivo points. 

10. An organism magnetic field-measured outcome display apparatus as set 
forth in claim 9, in which said reference graphic is in the form of a circular shape. 
[Detailed Explanation of the Invention] 

[0001] 

[Industrial Field in which the Invention is Used] 

The present invention relates to a tomogram image processing method, an 
organism magnetic field measuring method, an organism magnetic field-measured 
outcome display method, and apparatuses for carrying out the respective methods, and 
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more particularly, to a tomogram image processing method of and apparatus for 
obtaining fine image data based on tomogram images sliced at respective 
predetermined intervals, an organism magnetic field measuring method of and 
apparatus for analyzing an organism magnetic field source based on fine image data 
provided by the tomogram image processing method, and an organism magnetic 
field-measured outcome display method of and apparatus for displaying the analyzed 
outcome of the organism magnetic field source in an easily recognizable condition. 
[0002] 

[Description of the Prior Art] 

There have been provided a wide variety of conventional methods for visually 
superimposing magnetic field sources on tomogram images such as for example MRI 
images in order to display an organism magnetic field-measured outcome in an easily 
recognizable condition. As a method of visually superimposing magnetic field 
sources on tomogram images have been proposed a method of two-dimensionally 
displaying the magnetic field source under the condition that the magnetic field source 
is superimposed on the closest tomogram image and a method of three-dimensionally 
displaying the magnetic field source under the condition that the magnetic field source 
is superimposed on a plurality of tomogram images. 
[0003] 

As methods of measuring an organism magnetic field there have been 
proposed two methods consisting of a first measuring method and a second measuring 
method. The first measuring method is to be carried out under an assumption that the 
magnetic field source has a single dipole, and the second measuring method is to be 
carried out under an assumption that the magnetic field source has a plurality of current 
elements. The first measuring method comprises the step of calculating the position 
and the direction of the single dipole based on a plurality of outcomes of magnetic 
fields measurement measured at a plurality of points in the vicinity of a body surface. 
The second measurement method comprises the steps of measuring a magnetic field 
assuming that the magnetic field source has a predetermined number of current 
elements (usually 10 to 100); detecting points at which the measured magnetic field is 
remarkably greater than the predetermined number of current elements; and calculating 
the position and the direction of the current elements by performing a computation at 
the detected points. 
[0004] 

Furthermore, as methods of displaying an outcome of a measured organism 
magnetic field, there are publicly known an arrow mapping display method and a 
contour mapping display method. The arrow mapping display method comprises the 
steps of measuring an organism magnetic field, analyzing a magnetic field based on the 
measured organism magnetic field to obtain a plurality of current elements, and 
displaying the current elements thus obtained in the form of a current vector. The 
contour mapping display method comprises the steps of displaying the measured 
organism magnetic field in the form of an isomagnetic line. Furthermore, in order to 
carry out such display methods in a medical diagnoses, it is also proposed to display 
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arrows indicative of the current vectors or the isomagnetic lines superimposed on the 
MRI pictorial images so as to visually represent the current elements or the 
isomagnetic lines in internal organs. 
[0005] 

[Problem(s) to be Solved by the Invention] 

In the case that the two-dimensional display method in the abovementioned 
magnetic field source display methods is adopted, there is an inconvenience that it 
becomes very difficult for a medical doctor to recognize the position of the magnetic 
field source even though the medical doctor would carefully observe a plurality of the 
tomogram images by the reason that the tomogram image having the magnetic field 
source superimposed thereon is scarcely displayed. Moreover, although the 
tomogram image having the magnetic field source superimposed thereon should be 
selected and displayed, it is quite difficult to detect the relative position of the magnetic 
field source based on the selected tomogram image. 
[0006] 

Further, in the case that the three-dimensional display method in the 
abovementioned magnetic field source display methods is adopted, there is another 
inconvenience that it is very difficult to recognize how the magnetic field source is 
inwardly deep and widespread by the reason that there is no marks indicating the 
magnetic field except for the tomogram images. Moreover, it is also difficult to 
detect the position of the magnetic field relative to the internal organs because there are 
provided only a few tomogram images having the magnetic field superimposed thereon. 
In the case that the aforementioned first measuring method to be carried out under an 
assumption that the magnetic field source has a single dipole is adopted, there is 
further another inconvenience that the measured outcome is apt to converge on an 
excessively profound solution in addition to being not capable of measuring a 
distributive magnetic field source to be represented by a volume current flowing in the 
internal organ in pulsation by the reason that a single dipole is only obtained as the 
measured outcome. 
[0007] 

In the case that the aforementioned second measuring method to be carried out 
under an assumption that the magnetic field source includes a plurality of current 
elements is adopted, it becomes possible to measure the distributive magnetic field 
source. However, the fact that around 30-100 elements of currents are generally 
provided leads to the fact that it is required to calculate an estimate of each of current 
elements based on a plurality of magnetic field measuring values after setting up 
around 200-1000 channels of an array sensor and obtaining the plurality of magnetic 
field measuring values in the case that all of the current elements are to be measured. 
There is an inconvenience that the processing load to calculate the estimate in the 
abovementioned case requires T. FLOPS so that the system itself becomes 
considerably large in scale and cost. There is another inconvenience that the estimate 
thus calculated converges on the local minima, thereby remarkably aggravating the 
accuracy of the estimate of each of the current elements. 
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[0008] 

Any of the abovementioned arrow mapping display method and the contour 
mapping display method is appropriate for visually displaying the organism magnetic 
field-measured value in a particular timing, in other words, for visually displaying an 
organism magnetic field-measured distribution. This means that the both of the 
abovementioned arrow mapping display method and the contour mapping display 
method are inappropriate for displaying the organism magnetic field-measuring value 
of internal organs in pulsation. More particularly, although it is possible to 
dynamically display current elements by means of the arrow map display method, a 
large number of the current elements in the internal organs must be set up in order to 
display the volume current. Thus, a remarkably large number of the current elements 
must be displayed. Further, there is an inconvenience that it is difficult to recognize, 
for example, the propagation status of a visceral excitation current even though an MRI 
pictorial image and an arrow in an overlap display are carefully observed since the 
direction and dimension of each of the current elements change from hour to hour. 
[0009] 

The contour mapping display method cannot be adopted for continuously 
displaying the organism magnetic field-measured value by the reason that the 
isomagnetic line is unsuitable for displaying each of the current elements. This leads 
to the fact that the contour map display method cannot be carried out to visualize the 
propagation status of the visceral excitation current. 
[0010] 

[Purpose of the Invention] 

In order to solve the abovementioned problems, it is therefore a first object of 
the present invention to provide a tomogram image processing method and apparatus 
capable of obtaining fine image data based on fine tomogram image data of such as 
internal organs sliced at respective predetermined intervals, it is therefore a second 
object of the present invention to provide an organism magnetic field measuring 
method and apparatus capable of analyzing an organism magnetic field source easily 
and precisely based on the fine tomogram image data obtained by the tomogram image 
processing method, and it is therefore a third object of the present invention to provide 
an organism physical quantity-measured outcome display method and apparatus 
capable of continuously displaying the organism physical quantity-measured value in 
an easily recognizable condition. 
[0011] 

[Measures to Solve the Problems] 

In order to achieve the abovementioned first object, there is provided a 
tomogram image processing method as set forth in claim 1, comprising the steps of: a) 
carrying out an interpolation operation on a plurality of tomogram images of a test boy 
sliced in an area at predetermined intervals to obtain stereoscopic image data elements; 
b) dividing said area into small stereoscopic segments; c) counting the number of said 
stereoscopic image data elements for each of said small stereoscopic segments; and d) 
binarizing said small stereoscopic segments on the basis of the number of said 
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stereoscopic image data elements. 
[0012] 

A tomogram image processing apparatus as set forth in claim 2, comprises: 
interpolation carrying out means (3) for carrying out an interpolation operation on a 
plurality of tomogram images of a test boy sliced in an area at predetermined intervals 
to obtain stereoscopic image data elements; area dividing means (7) for dividing said 
area into small stereoscopic segments; stereoscopic image counting means (8) for 
counting the number of said stereoscopic image data elements for each of said small 
stereoscopic segments; and binarizing means (9, 10) for binarizing said small 
stereoscopic segments on the basis of the number of said stereoscopic image data 
elements. 
[0013] 

In order to achieve the foregoing second object, there is provided an organism 
magnetic field measuring method as set forth in claim 3, comprising the steps of: e) 
carrying out a plurality of Biot-Savart law calculations on a plurality of magnetic field 
information elements respectively for a plurality of vertexes of said small stereoscopic 
segments obtained in the tomogram image processing method as set forth in claim 1 ; f) 
adding up the results of the Biot-Savart law calculations to obtain a calculated 
magnetic field value; g) calculating a difference between said calculated magnetic field 
value and a measured magnetic field value; h) correcting data on said magnetic field 
sources included in the results of the Biot-Savart law calculations in accordance with 
said difference; i) judging whether or not said difference is vanishingly small; j) 
repeating said step (e) through step (h) when it is judged in said step (i) that said 
difference is not vanishingly small; and k) outputting, as an organism magnetic 
field-measured outcome, said data on said magnetic field sources included in the 
results of the Biot-Savart law calculations carried out in said step (h) when it is judged 
in said step (i) that said difference is vanishingly small. 
[0014] 

In an organism magnetic field measuring method as set forth in claim 4, said 
step (h) is replaced with the steps of: hi) calculating a difference between a coordinate 
of each of said vertexes of said small stereoscopic segments and a coordinate of each 
of corrected vertexes of said small stereoscopic segments; h2) calculating said 
coordinate of each of said corrected vertexes of said small stereoscopic segments on 
the basis of a code of said difference calculated in said step (hi); and h3) redistributing 
an estimate value of said magnetic field source in accordance with said difference 
calculated in said step (hi) for each of said corrected vertexes of said small 
stereoscopic segments calculated in said step (h2). 
[0015] 

An organism magnetic field measuring apparatus as set forth in claim 5, 
comprises: calculating means (111, 112, 11m; 101, 102, 10m) for carrying out 
a plurality of Biot-Savart law calculations on a plurality of magnetic field information 
elements respectively for a plurality of vertexes of said small stereoscopic segments 
obtained in the tomogram image processing method as set forth in claim 1; adding 
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means (20; 200) for adding up the results of the Biot-Savart law calculations outputted 
from said calculating means (111, 112, 11m; 101, 102, 10m) to obtain a 
calculated magnetic field value; error calculating means (30; 300) for calculating a 
difference between said calculated magnetic field value outputted from said adding 
means (20; 200) and a measured magnetic field value; correcting means (111a, 112a, 
lima) for correcting data on said magnetic field sources included the result of the 
Biot-Savart law calculations outputted from said calculating means (111, 112, 
11m; 101, 102, 10m) in accordance with said difference; and corrected result 
collecting means (40; 500) for collecting the results of the Biot-Savart law calculations 
corrected by said correcting means (111a, 112a, lima) and outputting them 
therethrough as an organism magnetic field-measured outcome. 
[0016] 

An organism magnetic field measuring apparatus as set forth in claim 6, 
further comprises: ,as correcting means, difference calculating means (41i) for 
calculating a difference between a coordinate of each of said vertexes of said small 
stereoscopic segments and a coordinate of each of corrected vertexes of said small 
stereoscopic segments; corrected vertex extracting means (42i) for calculating said 
coordinate of each of said corrected vertexes of said small stereoscopic segments on 
the basis of a code of said difference calculated by said difference calculating means 
(41i); and redistributing means (43i, 44i, 45i) for redistributing an estimate value of 
said magnetic field source in accordance with said difference calculated by said 
difference calculating means (41i) for each of said corrected vertexes of said small 
stereoscopic segments calculated by said corrected vertex extracting means (42i). 
[0017] 

In order to achieve the foregoing third object, there is provided a method as 
set forth in claim 7 of displaying an organism magnetic field-measured outcome along 
with an organism structure reference image, comprising the steps of: 1) calculating 
absolute values of said organism magnetic field-measured outcomes for a plurality of 
vivo points; m) setting up a size of a reference graphic for each of said vivo points on 
the basis of said absolute values calculated in said step (1); and n) displaying said 
reference graphic whose sizes is set up in said step (m) along with said organism 
structure reference image for each of said vivo points. In a method of displaying an 
organism magnetic field-measured outcome as set forth in claim 8, said reference 
graphic is in the form of a circular shape. 
[0018] 

An organism magnetic field-measured outcome display apparatus as set forth 
in claim 9 for displaying organism magnetic field-measured outcomes along with an 
organism structure reference image, comprises: reference graphic size calculating 
means (700) for calculating absolute values of said organism magnetic field-measured 
outcomes for a plurality of vivo points, and setting up a size of a reference graphic for 
each of said vivo points on the basis of said absolute values thus calculated; organism 
structure displaying means (800) for displaying said organism structure reference 
image; and reference graphic displaying means (900) for displaying said reference 
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graphic whose size is set up by said reference graphic size calculating means (700) for 

each of said vivo points. In an organism magnetic field-measured outcome display 

apparatus as set forth in claim 10, said reference graphic is in the form of a circular 

shape. 

[0019] 

[Function] 

A tomogram image processing method as set forth in claim 1, comprising the 
steps of: a) carrying out an interpolation operation on a plurality of tomogram images 
of a test boy sliced in an area at predetermined intervals to obtain stereoscopic image 
data elements; b) dividing said area into small stereoscopic segments; c) counting the 
number of said stereoscopic image data elements for each of said small stereoscopic 
segments; and d) binarizing said small stereoscopic segments on the basis of the 
number of said stereoscopic image data elements, can obtain the stereoscopic image 
data elements corresponding to all of the small stereoscopic segments by carrying out 
the interpolation operation on a plurality of tomogram images of the test body even 
through some tomogram images sliced at the intervals should fail to be provided. 
Further, the tomogram image processing method thus constructed can obtain a 
simplified stereoscopic image by the reason that the small stereoscopic segments are 
binarized on the basis of the number of the stereoscopic image data elements. This 
leads to the fact that the stereoscopic image obtained and displayed by the tomogram 
image processing method as set forth in claim 1 makes it possible for an operator to 
easily distinguish an area and/or a stereoscopic segment where an internal organ or the 
like is placed. 
[0020] 

Any of a linear interpolation arithmetical calculation, a spline interpolation 
arithmetical calculation, and the like may be appropriately selected and adopted for the 
abovementioned interpolation arithmetical operation. Further, an appropriate 
threshold value may be set up in order to adapt the numbers of the stereoscopic image 
data elements to be binarized to the numbers of the variety of the tomogram images 
and to the numbers of the internal organs to be processed. As a small stereoscopic 
segment, it is desirable to use a solid segment with which the segment can be filled 
densely. It more preferable to use a cube segment in order to simplify the treatments 
later performed. 
[0021] 

A tomogram image processing apparatus as set forth in claim 2, comprises: 
interpolation carrying out means (3) for carrying out an interpolation operation on a 
plurality of tomogram images of a test boy sliced in an area at predetermined intervals 
to obtain stereoscopic image data elements; area dividing means (7) for dividing said 
area into small stereoscopic segments; stereoscopic image counting means (8) for 
counting the number of said stereoscopic image data elements for each of said small 
stereoscopic segments; and binarizing means (9, 10) for binarizing said small 
stereoscopic segments on the basis of the number of said stereoscopic image data 
elements. The stereoscopic image obtained and displayed by the tomogram image 
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processing apparatus thus constructed makes it possible for an operator to easily 
distinguish an area and/or a stereoscopic segment where an internal organ or the like is 
placed. The organism magnetic field measuring method as set forth in claim 3, 
comprising the steps of: e) carrying out a plurality of Biot-Savart law calculations on a 
plurality of magnetic field information elements respectively for a plurality of vertexes 
of said small stereoscopic segments obtained in the tomogram image processing 
method as set forth in claim 1; f) adding up the results of the Biot-Savart law 
calculations to obtain a calculated magnetic field value; g) calculating a difference 
between said calculated magnetic field value and a measured magnetic field value; h) 
correcting data on said magnetic field sources included in the results of the Biot-Savart 
law calculations in accordance with said difference; i) judging whether or not said 
difference is vanishingly small; j) repeating said step (e) through step (h) when it is 
judged in said step (i) that said difference is not vanishingly small; and k) outputting, 
as an organism magnetic field-measured outcome, said data on said magnetic field 
sources included in the results of the Biot-Savart law calculations carried out in said 
step (h) when it is judged in said step (i) that said difference is vanishingly small, can 
obtain the data on the magnetic field source for each of a plurality of vertexes of said 
small stereoscopic segments by the reason that data on said magnetic field sources 
included in the results of the Biot-Savart law calculations is corrected until the 
difference between said calculated magnetic field value and a measured magnetic field 
value becomes vanishingly small. Furthermore, the number of unknown factors to be 
estimated can be largely reduced by the reason that the position of each of the vertexes 
of the small stereoscopic segments is known. 
[0022] 

The organism magnetic field measuring method as set forth in claim 4, in 
which said step (h) is replaced with the steps of: hi) calculating a difference between a 
coordinate of each of said vertexes of said small stereoscopic segments and a 
coordinate of each of corrected vertexes of said small stereoscopic segments; h2) 
calculating said coordinate of each of said corrected vertexes of said small stereoscopic 
segments on the basis of a code of said difference calculated in said step (hi); and h3) 
redistributing an estimate value of said magnetic field source in accordance with said 
difference calculated in said step (hi) for each of said corrected vertexes of said small 
stereoscopic segments calculated in said step (h2), can obtain accurate information on 
the magnetic field source by the reason that even magnetic field source vectors 
assigned in positions adjoining to each other, equal to each other in magnitude but 
opposing to each other in direction can be calibrated. 
[0023] 

In the organism magnetic field measuring apparatus as set forth in claim 5, the 
calculating means is operative to carry out a plurality of Biot-Savart law calculations 
on a plurality of magnetic field information elements respectively for a plurality of 
vertexes of said small stereoscopic segments obtained in the tomogram image 
processing method as set forth in claim 1, and the adding means is operative to add up 
the results of the Biot-Savart law calculations outputted from said calculating means to 
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obtain a calculated magnetic field value. The error calculating means is operative to 
calculate a difference between said calculated magnetic field value outputted from said 
adding means and a measured magnetic field value, and the correcting means is 
operative to correct data on said magnetic field sources included the result of the 
Biot-Savart law calculations outputted from said calculating means in accordance with 
said difference. After the correcting means has corrected the data on the magnetic 
field sources, the corrected result collecting means is operative to collect the results of 
the Biot-Savart law calculations thus corrected output them therethrough as an 
organism magnetic field-measured outcome. The organism magnetic field measuring 
apparatus thus constructed makes it possible for an operator to obtain the data on the 
magnetic field source for each of the vertexes of the small stereoscopic segments. 
Furthermore, the number of unknown factors to be estimated can be largely reduced by 
the reason that the position of each of the vertexes of the small stereoscopic segments 
is known. 
[0024] 

An organism magnetic field measuring apparatus as set forth in claim 6, 
further comprises: ,as correcting means, difference calculating means (41i) for 
calculating a difference between a coordinate of each of said vertexes of said small 
stereoscopic segments and a coordinate of each of corrected vertexes of said small 
stereoscopic segments; corrected vertex extracting means (42i) for calculating said 
coordinate of each of said corrected vertexes of said small stereoscopic segments on 
the basis of a code of said difference calculated by said difference calculating means 
(41i); and redistributing means (43i, 44i, 45i) for redistributing an estimate value of 
said magnetic field source in accordance with said difference calculated by said 
difference calculating means (41i) for each of said corrected vertexes of said small 
stereoscopic segments calculated by said corrected vertex extracting means (42i), can 
obtain accurate information on the magnetic field source by the reason that even 
magnetic field source vectors assigned in positions adjoining to each other, equal to 
each other in magnitude but opposing to each other in direction can be calibrated. 
[0025] 

The organism magnetic field measuring method will be described in detail 
hereinlater. As a result of enormous researches made by the inventor of the present 
invention in the field of measuring the organism magnetic field, it is found out that the 
position of each of the vertexes of the small stereoscopic segments is known because 
of the fact that each of the vertexes of the small stereoscopic segments has been 
obtained by the tomogram image processing method. This leads to the fact that 
information required to obtain the information on the magnetic field source is a vector 
component of the current element only. 
[0026] 

Even if the number of the small stereoscopic segments increases, the number 
of the vertexes does not increase in proportion to the increased number of the small 
stereoscopic segments by the reason that one vertex of one small stereoscopic segment 
serves as the vertex of the adjacent small stereoscopic segment. In the case that the 
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number of the small stereoscopic segments is large enough, the number of the vertexes 
and the number of the small stereoscopic segments are approximately equal to each 
other. By the fact that a current element of one vertex is related to that of an adjacent 
vertex, the amount of data of the current element of one vertex can be reduced 
remarkably, compared to the case that a general unknown quantity is estimated. 
[0027] 

The number of current elements to be finally estimated is remarkably small by 
the reason that the current element of a vertex having no magnetic field source 
converges on a zero vector. The inventions defined in claims 3 through 6 are made 
based on the abovementioned knowledge. A predetermined number of the magnetic 
field measuring values are obtained based on the efficiency and the capacity of the 
measuring system without considering an unknown number of current elements when 
the unknown current elements [the number is unknown either] are estimated at the start 
of measuring. A precise estimate of the unknown current elements can be achieved 
by repeating the correcting process of the data of the magnetic field source based on 
the predetermined number of the magnetic field measuring values, and the estimate 
outcome of which can be adopted as an organism magnetic field measured outcome. 
[0028] 

The organism magnetic field-measured outcome display method as set forth in 
claim 7 of displaying an organism magnetic field-measured outcome along with an 
organism structure reference image, comprising the steps of: 1) calculating absolute 
values of said organism magnetic field-measured outcomes for a plurality of vivo 
points; m) setting up a size of a reference graphic for each of said vivo points on the 
basis of said absolute values calculated in said step (1); and n) displaying said reference 
graphic whose sizes is set up in said step (m) along with said organism structure 
reference image for each of said vivo points, makes it possible for an operator to 
recognize, for example, the distributed condition of the maximal organism magnetic 
field-measured outcome in comparison with the case that the magnetic field measured 
outcome is displayed with the direction and the length of the arrows. The direction 
components of the organism magnetic field cannot be recognized. It is, however, not 
inconvenient by the reason that the organism magnetic field measured outcome can be 
dynamically displayed, thereby making it possible for an operator to easily recognize, 
for example, the time-variant fluctuation of the distributed condition of the maximal 
organism magnetic field measured outcome. The organism magnetic field-measured 
outcome display method thus constructed can satisfy the requirement of the medical 
analysis that the distribution of the maximal organism magnetic field-measured 
outcome and the time-variant fluctuation should be represented in an easily 
understandable manner. This leads to the fact that the organism magnetic 
field-measured outcome display method thus constructed is most appropriate for the 
medical analysis. 
[0029] 

In the organism magnetic field-measured outcome display method as set forth 
in claim 8, said reference graphic is in the form of a circular shape. In addition to the 
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operations as set forth in claim 7 the magnetic field-measured outcome display method 
as set forth in claim 8 can considerably enlarge a clearance between the reference 
graphics, thereby making it possible for an operator to recognize the organism structure 
reference image as a background. The organism magnetic field-measured outcome 
display apparatus as set forth in claim 9 is operative to display organism magnetic 
field-measured outcomes along with an organism structure reference image in a 
manner that the reference graphic size calculating means is operative to calculate the 
absolute values of said organism magnetic field-measured outcomes for a plurality of 
vivo points, and setting up a size of a reference graphic for each of said vivo points on 
the basis of said absolute values thus calculated, the organism structure displaying 
means is operative to display said organism structure reference image; and the 
reference graphic displaying means is operative to display said reference graphic 
whose size is set up by said reference graphic size calculating means for each of said 
vivo points. The organism magnetic field-measured outcome display apparatus thus 
constructed makes it possible for an operator to recognize, for example, the distributed 
condition of the maximal organism magnetic field-measured outcome in comparison 
with the case that the magnetic field measured outcome is displayed with the direction 
and the length of the arrows. The direction components of the organism magnetic 
field cannot be recognized. It is, however, not inconvenient by the reason that the 
organism magnetic field measured outcome can be dynamically displayed, thereby 
making it possible for an operator to easily recognize, for example, the time-variant 
fluctuation of the distributed condition of the maximal organism magnetic field 
measured outcome. The organism magnetic field-measured outcome display 
apparatus thus constructed can satisfy the requirement of the medical analysis that the 
distribution of the maximal organism magnetic field-measured outcome and the 
time-variant fluctuation should be represented in an easily understandable manner. 
This leads to the fact that the organism magnetic field-measured outcome display 
apparatus thus constructed is most appropriate for the medical analysis. 
[0030] 

In the organism magnetic field-measured outcome display apparatus as set 
forth in claim 10, said reference graphic is in the form of a circular shape. In addition 
to the operations as set forth in claim 9 the magnetic field-measured outcome display 
apparatus as set forth in claim 10 can considerably enlarge a clearance between 
reference graphics, thereby making it possible for an operator to recognize the 
organism structure reference image as a background. 
[0031] 

[EMBODIMENT OF THE PRESENT INVENTION] 

Preferred embodiments of the tomogram image processing method, the 
organism magnetic field measuring method, the organism magnetic field-measured 
outcome display method, and the apparatuses for the respective methods according to 
the present invention will now be described in detail in accordance with accompanying 
drawings. Referring now to the drawings, in particular to FIG. 1 of a flowchart, there 
is shown a preferred embodiment of the tomogram image processing method according 
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to the present invention. The tomogram image processing method comprises the 
steps of: waiting for receiving a plurality of tomogram images of a test boy sliced in an 
area at predetermined intervals (step SP 1); carrying out an interpolation operation on 
the plurality of received tomogram images to obtain interpolated tomogram images 
(step SP2); quantizing the received tomogram images and the interpolated tomogram 
images to obtain a plurality of squares each having a predetermined size (step SP3); 
dividing the area into a plurality of small stereoscopic segments each in the form of a 
cube shape having a center axis passing through the center of and a side equal in length 
to the side of the square (step SP4); counting the number of pixels constituting the 
received tomogram images and the interpolated tomogram images for each of the small 
stereoscopic segments (step SP5); judging whether or not the counted number is larger 
than a predetermined threshold value (for example, the average of all of the averages) 
for each of the small stereoscopic segments (step SP6), assigning each of the small 
stereoscopic segments, for which the counted number is judged in the step of SP6 to be 
larger than the threshold value, to "a small stereoscopic segment group to be 
displayed" (step SP7), and assigning each of the small stereoscopic segments , for 
which the counted number is judged in the step of SP7 to be not larger than the 
threshold value, to "a small stereoscopic segment group to be not displayed" (step SP8). 
The tomogram image processing method further comprises the steps of judging 
whether or not the step SP7 or the step SP8 has been performed for all of the small 
stereoscopic segments (step SP9). When it is judged in the step of SP9 that the step 
SP7 or the step SP8 has not been performed for a small stereoscopic segment, the step 
SP9 goes forward to the step SP10, in which the small stereoscopic segment is selected, 
and the step SP10 goes back to the step SP5. When it is judged in the step of SP9 that 
the step SP7 or the step SP8 has been performed for all of the small stereoscopic 
segment, the step SP9 goes forward to the step SP11, in which each of the small 
stereoscopic segments assigned to the small stereoscopic segment group to be 
displayed is visually displayed. 
[0032] 

It is needless to mention to the ratio of the number of pixels constituting the 
tomogram images to the number of pixels contained in a small stereoscopic segment 
for each of the small stereoscopic segments and a corresponding threshold value can be 
adopted in place of the number of pixels contained for each of the small stereoscopic 
segments. The present embodiment of the tomogram image processing method will 
be described in detail with reference to the drawings, in particular to FIGS. 2 through 4. 
FIG. 2 is a schematic diagram showing an example of a tomogram image, wherein the 
parts that are painted out indicates an internal organ of a human body. The tomogram 
image shown in FIG. 2 is constituted by a plurality of squares. As shown in FIG. 2, 
the center of each of the squares is shown with a small dot. This leads to the fact that 
a plurality of small stereoscopic segments each in the form of a cube shape having a 
center axis passing through the center of the square shown by the dot can be 
determined. As shown in FIG. 3, the small stereoscopic segments each having the 
number of pixels larger than a predetermined threshold value are assigned to the "small 
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stereoscopic segment group to be displayed" and the small stereoscopic segments each 
having the number of pixels not larger than a predetermined threshold value are 
assigned to the "small stereoscopic segment group to be not displayed". FIG. 4 is a 
perspective view corresponding to FIG. 3, and only the small stereoscopic segments 
assigned to the "small stereoscopic segment group to be displayed" is displayed. As 
will be seen in FIG. 4, the displayed small stereoscopic segments make it possible for 
an operator to visually recognize an internal organ. 
[0033] 

As will be understood from the foregoing description, carrying out a series of 
above-mentioned processes on all of the received tomogram images and the 
interpolated tomogram images makes it possible to visually display the organism 
internal organs as a solid model. It is, however, desirable to semi-transparently 
display the small stereoscopic segments assigned to the "small stereoscopic segment 
group to be displayed" so that the magnetic field source analyzing outcome can be 
easily recognized when a magnetic field source analyzing outcome is superimposed 
and displayed. Furthermore, specifying respective internal organs contained in the 
tomogram images makes it possible to display only the internal organs to be inspected 
or to display a plurality of the internal organs with their own colors. When an 
organism magnetic field source or the like, which will be described later, is analyzed, 
the tomogram image processing method according to the present invention thus 
constructed can process only the small stereoscopic segments assigned to the "small 
stereoscopic segment group to be displayed", thereby largely reducing arithmetical 
calculations required for the post processing by the reason that the small stereoscopic 
segments assigned to the "small stereoscopic segment group to be not displayed" 
correspond to a visceral cavity or the like and are accordingly ignored. 
[0034] 

[Second Embodiment] 

Referring to the drawings, in particular to FIG. 5 of a block diagram, there is 
shown a second embodiment of the tomogram image processing apparatus according 
to the present invention. The tomogram image processing apparatus, comprises: a 
tomogram image storing unit 1 for storing therein a plurality of tomogram images, a 
first tomogram image selecting unit 2 for selecting a tomogram image to be processed 
among the tomogram images stored in the tomogram image storing unit 1, an 
interpolation carrying out unit 3 for carrying out an interpolation operation on the basis 
of the tomogram images selected by the first tomogram image selecting unit 2 to obtain 
interpolated tomogram images, an interpolation tomogram image storing unit 4 for 
storing the interpolated tomogram images obtained by the interpolation carrying out 
unit 3, a first repeat control unit 5 for having the tomogram image selecting unit 2 and 
the tomogram image interpolating unit 3 repeatedly operate for a necessary number of 
times, a second tomogram image selecting unit 6 for sequentially selecting a tomogram 
image from among the tomogram images stored in the tomogram image storing unit 1 
or the interpolation tomogram image storing unit 4 in response to an end signal 
indicative of the end of the repeat control operation of the first repeat control unit 5, a 
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stereoscopic segment assigning unit 7 for assigning the tomogram image selected by 
the second tomogram image selecting unit 6 to a stereoscopic segment, a pixel number 
counting unit 8 for counting the number of pixels present in the stereoscopic segment, 
a judging unit 9 for judging whether or not the number of pixels counted by the pixel 
number counting unit 8 is larger than a predetermined threshold value, a 
display-stereoscopic segment assigning unit 10 for assigning the small stereoscopic 
segment, for which the counted number is judged by the judging unit 8 to be larger 
than the threshold value, to "a small stereoscopic segment group to be displayed", a 
non-display-stereoscopic segment assigning unit 11 for assigning the small 
stereoscopic segment, for which the counted number is judged by the judging unit 8 to 
be not larger than the threshold value, to "a small stereoscopic segment group to be not 
displayed", an assigning outcome storing unit 12 for storing therein the received 
tomogram images and the interpolated tomogram images in association with the small 
stereoscopic segments assigned to the respective small stereoscopic segment groups, a 
second repeat control unit 13 for controlling the second tomogram image selecting unit 
6, the stereoscopic segment assigning unit 7, the pixel counting unit 8, the judging unit 
9, the display-stereoscopic segment assigning unit 10 and the non-displaying assigning 
unit 11 to have the second tomogram image selecting unit 6, the stereoscopic segment 
assigning unit 7, the pixel counting unit 8, the judging unit 9, the display-stereoscopic 
segment assigning unit 10 and the non-displaying assigning unit 11 repeatedly operate 
for a necessary number of times, and a stereoscopic segment displaying unit 14 for 
visually displaying only the small stereoscopic segments assigned to the "small 
stereoscopic segment group to be displayed" in response to an end signal indicative of 
the end of the repeat control operation of the second repeat control unit 13. 
[0035] 

As will be seen from the foregoing description, it is to be understood that the 
tomogram processing apparatus according to the present invention thus constructed as 
shown in FIG. 5 can perform processing on the tomogram images as shown in FIGS. 2 
through 4, and thus visually display only the small stereoscopic segments assigned to 
the "small stereoscopic segment group to be displayed". 
[0036] 

[Third Embodiment] 

Referring to the drawings, in particular to FIG. 6 of a block diagram, there is 
shown a third preferred embodiment of the organism magnetic field measuring 
apparatus according to the present invention. The organism magnetic field measuring 
apparatus is shown in FIG. 6 as comprising: a plurality of Biot-Savart arithmetical 
calculating units 111, 112, ... ,11m each for carrying out a Biot-Savart law calculation 
on a magnetic field information element for a vertex of the small stereoscopic segment 
obtained in the tomogram image processing method shown in FIG. 1 or the tomogram 
image processing apparatus shown in FIG. 5, an adding unit 20 for adding up the 
results of the Biot-Savart law calculations outputted from the Biot-Savart arithmetical 
calculating unit 111, 112, 11m; 101, 102, an error calculator 30 for calculating 
a difference between an added-up value Oj(t) outputted from the adding unit 20 and a 
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magnetic field-measured value Sj(t) inputted as a teacher pattern, a plurality of 
correcting units 11 la, 112a, lima respectively for correcting vector components of 
current elements estimated by the Biot-Savart arithmetical calculating units 111, 
112, 11m based on the difference calculated by the error calculator 30, and a data 
collecting unit 40 for collecting the vector components of the current elements 
estimated by the Biot-Savart arithmetical calculating units 111, 112, 11m and to 
output the vector components thus collected as an analyzing outcome. Furthermore, 
each of the Biot-Savart arithmetical calculating units 111, 112, ... ,11m is adapted to 
carry out a Biot-Savart law calculation on a magnetic field information element for a 
vertex of the small stereoscopic segment obtained at time "t" in the tomogram image 
processing method shown in FIG. 1 or by the tomogram image processing apparatus 
shown in FIG. 5, and correct the vector components of the estimated current elements 
in response to the difference dj(t) = {Sj(t) - Oj(t)} calculated by the error calculator 30. 
The operations of the Biot-Savart arithmetical calculating units 111, 112, ... ,11m may 
be performed synchronously or asynchronously. 
[0037] 

The operation of the organism magnetic field measuring apparatus thus 
constructed will be described hereinlater. Magnetic field Oj(t) to be analyzed is 
represented by the linear sum of the time "t", three-dimensional coordinate values xi, yi, 
zi (" i" is intended to mean a positive integer) of each vertex of the respective small 
stereoscopic segments, and functions gi (each of the functions is determined based on 
the Biot-Savart law). Each of the functions gi has vector components Pxi, Pyi, Pzi of 
the current elements. When only vector components on a plane are used as 
measurable components for a magnetic field sensor, any of the vector components can 
be omitted by applying the coordinates to the sensor. All of the three-dimensional 
coordinate values are known by the fact that the respective small stereoscopic 
segments are provided by the tomogram image processing method of FIG. 1 or the 
tomogram processing apparatus of FIG. 5. 
[0038] 

Accordingly, the organism magnetic field measuring apparatus thus 
constructed can obtain the added-up value Oj(t) with the processes of: supplying "m" 
unit of the Biot-Savart arithmetical calculating units 111, 112, ... , 11m with the time 
"t" and the three-dimensional coordinate values of each of the vertexes of the 
respective small stereoscopic segments; performing the arithmetical calculation gl, 
g2, gm respectively to calculate the function values, and supplying the adding unit 
20 with all of the calculated values. The added-up value O j (t) provided on the first 
iteration of the above-mentioned processes is however different from the actual 
measured value S j (t) by the reason that the vector components of the respective 
current elements are set up at arbitrary numbers on the first iteration. The error 
calculator 30 is then operated to calculate a difference between actual measured value 
Sj(t) and the added-up value O j (t), and feed the difference thus calculated as an 
estimated error dj (t) back to the respective correcting units 111a, 112a, ... , lima so 
that the Biot-Savart arithmetical calculating units 111, 112, ... , 11m will alter the 
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vector components of the respective Biot-Savart arithmetical calculating units until the 

estimated error dj (t) becomes small. 

[0039] 

The above-mentioned processes are repeated until the estimated error dj (t) 
becomes approximately zero. The data collecting unit 40 is then operated to collect 
the vector components of the current elements estimated by the Biot-Savart 
arithmetical calculating units 111, 112, 11m and to output the vector components 
thus collected as an analyzing outcome. The vector component of the current element 
in each of the vertexes of the respective small stereoscopic segments can be outputted 
from the data collecting unit 40 at this stage. Further, an estimated error evaluating 
function Ej (t) is defined by the formula: Ej(t) = (1/2) {Sj (t) - Oj (t) } 2 , Equation 1 is 
obtained as follows. 
[Equation 1] 

3Ej (t) /aoj (t) — [S j (t) -Oj (t) ) 

[0041] 

If the vector component estimated by each of the Biot-Savart arithmetical 
calculating units is corrected based on a steepest descent method, the vector 
component, which minimizes the value of the estimated error evaluating function, can 

be obtained in accordance with the following Equation 2. Here, e k is a learning gain 

(a correcting gain) of the vector component aik. 

[0042] 

[Equation 2] 

a i k = a i k - e k (9E j (t) /3a i kl 
■alk-fk l(9Ej (t) /90 j (t)} (90 j CO/ 
9aik] - a i k + e k {Sj (t)-Oj (O) { 9 O j CO 
/9a ik) 

[0043] 

The equation 2 can be expressed with the equation 4 by the reason that the 
partial differential value with aik of the added-up value O j (t) is given in accordance 
with the equation 3. 
[0044] 
[Equation 3] 

90 j (t)/3aik = 3 ||g i (t, ail, a i 2, a i 3, 
• • • a i L) ) / da i k 

- 9 IB i ( t. a i 1, a I 2, a i 3. •••aiL)) / 

9 a i k 

[0045] 
[Equation 4] 

aik = aikH-ek {Sj (t) -Oj (t)} [9 (g i ( t , 
ail, ai2, a i 3, • ■ ■ a i L) ] /9 a L k] 
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[0046] 

From the foregoing description, it is to be understood that the second 
preferred embodiment of the present invention, which carries out the operation in 
accordance with Equation 4, can thus improve the accuracy of the estimated vector 
components. In addition, the corrected value may be set at a negative value when the 
gradient of the estimated error dj(t) is positive, and may be set at a positive value when 
the gradient of the estimated error dj(t) is negative. When the vector component of 
the current element of each vertex is thus obtained, a vector space interpolation 
processing is carried out based on the current elements of eight vertexes, a magnetic 
flux density where a minute point occurs is obtained, the magnetic flux density where 
each of the minute points occurs is integrated, and a magnetic field of each of the small 
stereoscopic segments is then obtained. Thus, the distributed currents can be 
displayed with fogs or arrows based on the provided outcome. 
[0047] 

Furthermore, the process performed by the apparatus shown in FIG. 6 is 
applicable only to the small stereoscopic segments assigned to the "small stereoscopic 
segment group to be displayed" obtained by the embodiment shown in FIG. 1 or FIG. 
5, thereby largely reducing arithmetical calculations by the reason that the small 
stereoscopic segments assigned to the "small stereoscopic segment group to be not 
displayed" correspond to a visceral cavity or the like and are accordingly ignored. 
Moreover, even if the number of the small stereoscopic segments increases, the number 
of the vertexes does not increase in proportion to the increased number of the small 
stereoscopic segments by the reason that one vertex of one small stereoscopic segment 
serves as the vertex of the adjacent small stereoscopic segment. In the case that the 
number of the small stereoscopic segments is large enough, the number of the vertexes 
and the number of the small stereoscopic segments are approximately equal to each 
other. This also leads to the fact that the present embodiment can reduce the 
operation load. 
[0048] 

According to the present invention, the vector components can be sufficiently 
estimated with a comparatively small number of the measuring points although it has 
been generally believed that the number of measuring points is required to be greater 
than the number of unknown values in order to obtain unknown magnetic 
field-measured values. This also leads to the fact that the present embodiment can 
degrease the number of measuring points, thereby as a whole reducing the operation 
load. 
[0049] 

The following description will be directed to the above-mentioned aspect in 
detail. Conventionally, internal organs has not been defined as aggregations of small 
stereoscopic segments, and, for example, the number of measuring points is required to 
be not less than the number of the current elements to be estimated. This results in 
the fact that number of the measuring points tends to be more than required especially 
when the number of the current elements is unknown. The number of the current 
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elements could not but become enormous by the reason that the currents flowing in the 
internal organs vary in position, direction, and absolute value with the course of time. 
In the present embodiment, on the other hand, since the position of each of the small 
stereoscopic segments is fixed in advance; the current element for a small stereoscopic 
segment, which does not influence the measuring point, converges on a zero vector. 
The current element, which has converged on the zero vector, is not required to be 
processed later. This leads to the fact that the present embodiment of the organism 
magnetic field measuring apparatus can accurately estimate the vector components 
based on magnetic field-measured values obtained at a relatively small number of 
measuring points. 
[0050] 

Referring now to the drawings, in particular to FIG. 7 of a block diagram, 
there is shown another preferred embodiment of the organism magnetic field 
measuring apparatus according to the present invention, and in particular to FIG. 8 of a 
schematic diagram, there is shown a detailed construction of a part corresponding to a 
hierarchical type perceptron. The present embodiment of the organism magnetic field 
measuring apparatus is as shown in FIGS. 7 and 8 comprising: a plurality of 
hierarchical type perceptrons lOp (p=l,2, ... ,m) each having a plurality of inputs 
terminals and one output terminal, an adding unit 200 for adding up function values gij 
outputted from the respective hierarchical type perceptrons lOp, a function value error 
calculator 300 for calculating a difference between the added-up outcome Oj(t) 
outputted from the adding unit 200 and a measuring value Sj(t) inputted as a teacher 
pattern, a plurality of partial differential value calculating units lOpr (r=l,2, ... n) each 
for calculating a partial differential value corresponding to the function value gij 
outputted from the hierarchical type perceptron lOp, a partial differential value error 
calculator 30pr each for calculating a difference between an output from the partial 
differential value calculating unit lOpr of the hierarchical type perceptron lOp and a 
partial differential value (a partial differential value teacher pattern) calculated in 
advance in accordance with the numerical differentiation method, a partial differential 
function learning unit 60p each for having the applicable hierarchical type perceptron 
lOp learn a partial differential function based on the difference calculated by the 
function value error calculator 300 and the difference calculated by the partial 
differential value error calculator 30pr, a plurality of correcting units 12p each for 
correcting data to be inputted into the input layer of the hierarchical type perceptron 
lOp (hereinlater referred to as "input data of the hierarchical type perceptron lOp") 
based on the difference between the added-up outcome Oj(t) calculated by the function 
value error calculator 300 and the measured value Sj(t) as the teacher pattern, and 
based on the partial differential value calculated by the partial differential value 
calculating unit lOpr, a control unit 400 for selecting a correcting unit 12p and a partial 
differential function learning unit 60p from among the plurality of the correcting units 
12p and partial differential function learning units 60p and controlling the correcting 
unit 12p and the partial differential function learning unit 60p thus selected to have the 
selected correcting unit 12p and the selected partial differential function learning unit 
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60p operate repeatedly for a predetermined number of times (until the difference 
becomes small enough), and a data collecting unit 500 for outputting the input data of 
the respective hierarchical type perceptrons lOp corresponding to unknown vector 
components as unknown quantity estimated outcomes under the condition that the 
input data of the hierarchical type perceptrons lOp has been corrected by the correcting 
units 12p for the predetermined number of times. It is herein assumed that each of 
the hierarchical type perceptrons lOp has sufficiently learned functions based on the 
Biot-Savart law with known input patterns and teacher patterns. Furthermore, the 
operations of the hierarchical type perceptrons lOp may be performed synchronously 
or asynchronously. 
[0051] 

Furthermore, each of the partial differential function learning units 60p is 
adapted to alternately repeat the two ways of learning: namely, the learning of the 
hierarchical type perception lOp based on the difference calculated by the function 
value error calculator 300 (for example, a learning of back propagation) and the 
learning of the hierarchical type perception lOp based on the difference calculated by 
the partial differential value error calculator 30pr, thereby making it possible for each 
of the hierarchical type perceptrons lOp to learn the partial differential function 
without losing the learning results so much. 
[0052] 

Each of the correcting units 12p is adapted to correct an unknown quantity 
based on the difference between the added-up outcome Oj(t) calculated by the function 
value error calculator 300 and the measured value Sj(t) as the teacher pattern, and 
based on the partial differential value calculated by the partial differential value 
calculating unit lOpr. The calculating unit lOpr, correcting unit 12p, and the function 
value error calculator 300 collectively constitute an unknown quantity estimating 
section. The operation of the unknown quantity estimating section will be described 
in detail hereinlater. The unknown quantity estimating section is operated to slightly 
change the input data of one of the hierarchical type perceptrons lOp corresponding to 
an unknown vector component after the hierarchical type perceptrons lOp have 
completed the learning based on the Biot-Savart law in order to obtain an output 
pattern and a corresponding partial differential value. Also, it is assumed that the 
unknown quantity estimating apparatus has calculated in advance a partial differential 
value in the case that the known input corresponding to the unknown vector 
component has been slightly changed in accordance with the numerical differentiation 
method or the like, and stored therein the partial differential value as a partial 
differential value teacher pattern. In the condition thus prepared, the control unit 400 
is operated to select a partial differential function learning unit 60p from a plurality of 
partial differential function learning units 60p and control the selected partial 
differential function learning unit 60p to have the selected partial differential function 
learning unit 60p operate repeatedly for a predetermined number of times. In the way, 
a back propagation learning is carried out, and a load and a threshold value of a 
neuronal component constituting each of the hierarchical type perceptrons lOp is 
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determined so that a function based on the Biot-Savart law and an operation equivalent 
to the partial differential numerical arithmetic of the function can be realized. 
[0053] 

Furthermore, the load and the threshold value thus determined as a learning 
outcome for the selected hierarchical type perceptrons lOp can be adopted as a load 
and a threshold value for each of hierarchical type perceptrons other than the selected 
hierarchical type perceptrons lOp. This leads to the fact that the hierarchical type 
perceptrons other than the selected hierarchical type perceptrons lOp are not required 
to learn their partial differential functions, thereby drastically reducing the learning 
time. After the learning process is completed, in place of the partial differential 
function learning unit 60p, the control unit 400 is operated to select a correcting unit 
12p from among a plurality of correcting units. The data collecting unit 500 is then 
operated to output the input data of the respective hierarchical type perceptrons lOp 
corresponding to unknown vector components as unknown quantity estimated 
outcomes when the input data of the hierarchical type perceptrons lOp has been 
corrected by the correcting units 12p for the predetermined number of times. As will 
be seen from the foregoing description, it is to be understood that the present 
embodiment of the organism magnetic field measuring apparatus can estimate the 
vector components as accurately as the embodiment shown in FIG. 6. 
[0054] 

[Fourth Embodiment] 

Referring to the drawings, in particular to FIG. 9 of a block diagram, there is 
shown a fourth embodiment of the organism magnetic field measuring apparatus 
according to this present invention. The constitutional elements of the present 
embodiment of the organism magnetic field measuring apparatus as shown in FIG. 9 
are the entirely the same as those of the previous embodiment of the organism 
magnetic field measuring apparatus according to this present invention as shown in 
FIGS. 7 and 8 except for the fact that the present embodiment of the organism 
magnetic field measuring apparatus comprises an input changing unit 15p for slightly 
changing an input, an output storing unit 16p for storing therein an output in response 
to the input not changed by input changing unit 15p and an output in response to the 
input changed by the input changing unit 15p, and a difference calculating unit 17p for 
calculating the difference between the output in response to the input not changed by 
input changing unit 15p and the output in response to the input changed by the input 
changing unit 15p, respectively stored in the output storing unit 16p, in place of the 
partial differential value calculating units lOpr provided in the respective hierarchical 
type perceptrons in the previous embodiment shown in FIGS. 7 and 8. 
[0055] 

Accordingly, the present embodiment of the organism magnetic field 
measuring apparatus can calculate the difference between the output in response to the 
regular input and the output in response to the slightly changed input. The difference 
between the output in response to the regular input and the output in response to the 
slightly changed input corresponds to a partial differential value. This leads to the 
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fact that the present embodiment of the organism magnetic field measuring apparatus 
can achieve the same advantage of the previous embodiment shown in FIGS. 7 and 8 
while eliminating the need of learning the partial differential function, thereby largely 
reducing the leaning time. 
[0056] 

[Fifth Embodiment] 

Referring to the drawings, in particular to FIG. 10 of a block diagram, there is 
shown a unit of a fifth embodiment of the organism magnetic field measuring 
apparatus according to this present invention. The organism magnetic field 
measuring apparatus is shown in FIG. 10 as comprising: a plurality of Biot-Savart 
arithmetical calculating units Hi each for performing an arithmetical calculation in 
accordance with the Biot-Savart law, a current element collecting unit 40a for 
collecting vector components of current elements estimated by the Biot-Savart 
arithmetical calculating units Hi and for outputting the collected vector components as 
an analysis outcome, a plurality of correction quantity calculating units 41i each for 

calculating correction quantities Ax, Ay, and Az on the basis of the 

three-dimensional coordinate values provided by the method shown in FIG. 1 and/or 
by the apparatus shown in FIG. 5, a plurality of vertex extracting units 42i each for 
extracting coordinate of each of the adjacent vertexes, for which the estimated vector 
component of the current element is to be distributed, the basis of the code of the 

calculated correction quantities Ax, Ay, and Az, a plurality of standardizing unit 43i 

each for standardizing a correction vector calculated based on the correction quantities 

Ax, Ay, and Az, a plurality of redistributing ratio calculating unit 44i each for 

calculating a redistributing ratio based on the absolute value of each component of the 
standardized correction vectors, a plurality of redistributing units 45i each for 
redistributing an estimated value to the vector component of the current element 
estimated for each of the extracted vertexes on the basis of the estimated vector 
component of the current element estimated and the calculated redistributing ratio, a 
plurality of estimated value updating units 46i each for updating the estimated vector 
component of the current element in response to the redistributed estimated value, and 
a repeat control unit 47 for controlling the Biot-Savart arithmetical calculating unit Hi, 
the correction quantity calculating unit 41i, the vertex extracting unit 42i, the 
standardizing unit 43i, the redistributing ratio calculating unit 44i, the redistributing 
unit 45i and the estimated value updating unit 46i to have the Biot-Savart arithmetical 
calculating unit Hi, the correction quantity calculating unit 41i, the vertex extracting 
unit 42i, the standardizing unit 43i, the redistributing ratio calculating unit 44i, the 
redistributing unit 45i and the estimated value updating unit 46i repeatedly operate 

until the correction quantities Ax, Ay, and A z becomes vanishingly small. 

[0057] 



23 



The estimated value updating unit 46i is operative to update the estimated 
vector component of the current not by simply adding the redistribution ratio to the 
estimated vector component of the current element, but by adding the average of all of 
the redistribution ratios to the estimated vector component of the current element. 
The correction quantity calculating unit 41i is operative to calculate the correction 

quantities Ax, Ay, Az by computing the partial differential value of the function 

value calculated by for example Biot-Savart arithmetical calculation units. 



The operation of the organism magnetic field measuring apparatus thus 
constructed will be described hereinlater. Each of the Biot-Savart arithmetical 
calculation units Hi is operated to carry out the computation in accordance with the 
following equation: 



where i 0 is intended to mean a magnetic permeability, i (i =1,2, ... , N) is intended to 
mean the current element contained in each vertex of the respective small stereoscopic 
segments, j (j=l,2, ... ,m) is intended to mean measuring points. When it is assumed 
that the measuring condition of each of the measuring points is represented by Pj = (xj, 
yj, zj) and the unknown physical quantity of each of the current elements is represented 
by Ui = (xi, yi, zi, pxi, pyi), the following equation is derived. 

Xji = xj - xi, Yji = yj - yi, Zji = zj - zi, and Rji = (Xji 2 + Yji 2 + zji 2 ) m . 



Thus, the partial differential values of the above-mentioned equations partially 
differentiated with the respective unknown quantities can be obtained by supplying the 
correction quantity calculating units 41 with the outcome respectively outputted from 
the Biot-Savart arithmetical calculating units Hi. The partial differential values can 
be computed in accordance with the following Equations 5 through 9. 
[0060] 
[Equation 5] 



[0058] 



gij = (i 0/4 6) {(Yji x pxi -Xji x pyi) /Rji 3 } 



[0059] 



3 



8gij 




[0061] 
[Equation 6] 




[0062] 
[Equation 7] 




[0063] 
[Equation 8] 
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[0064] 
[Equation 9] 



9g i j 

9Py i 
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X j i 



R j i 3 



[0065] 



The correction quantities Ax, Ay, Az can be obtained based on the partial 

differentiation computation carried out in accordance with the Equations 5 through 7 

and Sj— Oj. With the correction quantities Ax, Ay, A z thus obtained, the adjacent 

vertexes, for which the estimated vector components of the current elements are to be 
distributed, can be extracted by the vertex extracting unit 42i based on the codes of the 

correction quantities Ax, Ay, Az. As best shown in FIG. 11, concretely, vertexes 

respectively shown by double circles are extracted as adjacent vertexes, for which the 
estimated vector components of the current elements are to be distributed, by the fact 
that the correction vectors respectively shown by dotted lines are obtained. The 
correction vectors are then standardized by the standardizing unit 43i, the redistributing 
ratio is calculated by the redistributing ratio calculating unit 44i based on the absolute 
value of each component of the standardized correction vectors, and an estimate is 
redistributed by the redistributing unit 45i to the vector component of the current 
element estimated for each of the extracted vertexes extracted by the extracting unit 
42i. In this case, each of the vector components of the current elements may be 
drastically redistributed on the basis of a plurality of correction vectors. Simply 
adding the vector components of the current elements may result in the fact that the 
value before the redistribution is by far different from the value after the redistribution, 
thereby undemining the reliability of the stability of the estimation. In order to solve 
this problem, the estimated value updating unit 46i is operative to update the estimated 
vector component of the current by adding to the estimated vector component of the 
current element the average of all of the redistribution ratios calculated based on a 
plurality of standardized correction vectors. 



The repeat control unit 47 is operated to control the Biot-Savart arithmetical 
calculating unit Hi, the correction quantity calculating unit 41i, the vertex extracting 
unit 42i, the standardizing unit 43i, the redistributing ratio calculating unit 44i, the 
redistributing unit 45i and the estimated value updating unit 46i to have the Biot-Savart 
arithmetical calculating unit Hi, the correction quantity calculating unit 41i, the vertex 
extracting unit 42i, the standardizing unit 43i, the redistributing ratio calculating unit 
44i, the redistributing unit 45i and the estimated value updating unit 46i repeatedly 

operate until the correction quantities Ax, Ay, and A z becomes vanishingly small for 



[0066] 
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all of the vertexes, thereby highly improve the accuracy of the estimates of the current 
elements. While the estimation process is performed by the repeat control unit 47, 
Equation 5, Equation 6, and Equation 7 are not directly influenced by the correction 

quantities Ax, Ay, A z by the reason that the three-dimensional coordinate values of 

each vertex of the respective small stereoscopic segments are to be fixed. The 
aforementioned equations are, on the other hand, indirectly influenced by the 

correction quantities Ax, Ay, A z by the fact that the vector components of the current 

elements are changed through the redistribution process. This means that through out 
the processes carried out in accordance with Equation 5, Equation 6, Equation 7, 
Equation 8, and Equation 9, Xji, Yji, Zji, and Rji, remain unchanged while, on the 
other hand, only pxi and pyi are changed, thereby ensuring that the highly estimated 
pxi and pyi are obtained. The conventional organism magnetic field measuring 
apparatus tends to remain currents canceling each other at positions at which in fact no 
current flows even if the currents canceling each other are placed in the vicinity of 
each other as shown, for example, in FIG. 12(A). The present embodiment of the 
organism magnetic field measuring apparatus, on the other hand, carries out the 
redistribution process, thereby making it possible for the currents canceling each other 
to be displayed at positions at which actually the currents flow as shown in FIG. 12(B). 
In FIGS. 12(A) and 12(B), a black dot is intended to mean a vertex to be processed and 
a circle is intended to mean a vertex placed in the vicinity of the vertex. 
[0067] 

Examples of the magnetic field source estimating outcomes obtained 
approximately at the time the Q wave occurs based on magnetic field distribution maps 
are shown in FIGS. 14 to 18. The magnetic field distribution maps have been 
prepared by synchronously adding a plurality of magnetocardiograms measured at 3 
cm intervals with 6x6 points to an electrocardiography using a one-channel magneto 
meter. As will be seen from FIGS. 14 to 18, it is to be understood that the magnetic 
field source estimating outcomes precisely represent a cardiac volume current, which is 
medically confirmed, with a high degree of accuracy. FIG. 13 show information on 
vertexes with depth of 6 cm obtained based on MM pictorial images. FIGS. 14 to 18 
show estimating outcomes respectively obtained at the time 35 msec, 33 msec, 30 
msec, 25 msec, and 20 msec prior to the occurrence of R wave. The large black dots 
collectively represent a cardiac wall, and line segments respectively originating from 
the black dots represent current elements. 
[0068] 

In the present embodiment of the organism magnetic field measuring 
apparatus according to the present invention, the hierarchical type perceptrons may be 
adopted in place of the Biot-Savart arithmetical calculating units. 
[0069] 

[Seventh Embodiment] 

Referring to the drawings, in particular to FIG. 19 of a flowchart, there is 
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shown a seventh preferred embodiment of the organism magnetic field measured 
outcome display method according to this present invention. As shown in FIG. 19, 
the organism magnetic field measured outcome display method comprises the steps of: 
waiting for receiving the position and the magnitude (the absolute value of an 
excitation current) of the current element in a plane corresponding to the tomogram 
image at a predetermined point of time (step SP 1), assigning a circle having a 
diameter corresponding to the magnitude of each of a plurality of the current elements 
to each of the current elements (step SP 2), displaying the organism tomogram image 
(step SP 3), displaying the circle assigned to each of the current elements along with 
the organism tomogram image for each of vivo points (step SP 4), judging whether or 
not the end of the excitation current display is instructed (step SP 5), and waiting for 
receiving the position and the magnitude of the current element in a plane 
corresponding to the tomogram image at another point of time when it is judged in the 
step SP5 that the end of the excitation current display is not instructed (step SP 6). 
Then, the step SP6 goes back to the step SP2. When it is, on the other hand, judged 
in the step SP5 that the end of the excitation current display is not instructed, the 
process goes to END. In the step SP2, the diameter of the circle assigned to each of 
the current elements may precisely correspond to the magnitude of each of the current 
elements. Alternatively, the magnitude of each of the current elements may be 
divided into a plurality of magnitude grades, and the diameter of the circle assigned to 
each of the current elements may correspond to the magnitude grade of each of the 
current elements. 
[0070] 

[Eighth Embodiment] 

Referring to the drawings, in particular to FIG. 20 of a flowchart, there is 
shown an eighth preferred embodiment of the organism magnetic field measured 
outcome display method according to this present invention. The organism magnetic 
field measured outcome display method comprises the steps of: waiting until a point at 
which an internal organ is placed is detected on a plane corresponding to the tomogram 
image (step SP 1), assigning a circle having a predetermined diameter to the detected 
point at which the internal organ is placed (step SP 2), waiting until the position and 
the magnitude of each of the current elements is obtained at a predetermined point of 
time on the above-mentioned plane (step SP 3), replacing the circle assigned to the 
point in the step SP2 with a high-lightened circle having a diameter corresponding to 
the magnitude of each of the current elements (step SP 4), displaying the organism 
tomogram image (step SP 5), displaying the circle assigned to the detected point along 
with the organism tomogram image at the detected point (step SP 6), judging whether 
or not the end of the excitation current display is instructed (step SP 7), and waiting for 
the position and the magnitude of the current element in a plane corresponding to the 
tomogram image at another point of time when it is judged in the step SP7 that the end 
of the excitation current display is not instructed (step SP 8). Then, the step SP8 goes 
back to the step SP4. When it is, on the other hand, judged in the step SP7 that the 
end of the excitation current display is not instructed, the process goes to END. 
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[0071] 

Examples of the excitation currents displayed by the present embodiment of 
the organism magnetic field measured outcome display method are shown in FIG. 21. 
As will be seen from FIG. 21, the excitation currents are remarkably visible. The 
circles are displaced at the point at which the internal organ is placed. In addition, a 
clearance between the circles is enlarged, thereby making it possible for an operator to 
recognize the tomogram image as a background. In particular, the present 
embodiment of the organism magnetic field measured outcome display method can 
display the excitation currents on the time series, thereby making it possible for an 
operator to easily recognize, for example, the condition of the maximal excitation 
current in the process of changing. 
[0072] 

[Ninth Embodiment] 

Referring to the drawings, in particular to FIG. 22 of a flowchart, there is 
shown a ninth preferred embodiment of the organism magnetic field measured 
outcome display method according to this present invention. The organism magnetic 
field measured outcome display method comprises the steps of: waiting for receiving 
the position and the magnitude of the current element at a predetermined point of time 
(step SP 1), assigning a sphere having a diameter corresponding to the magnitude of 
each of a plurality of the current elements to each of the current elements (step SP 2), 
three-dimensionally displaying a semi-transparent organism structure reference image 
(step SP 3), displaying the sphere assigned to each of the current elements along with 
the semi-transparent organism structure reference image for each of corresponding 
points (step SP 4), judging whether or not the end of the excitation current display is 
instructed (step SP 5), and waiting for receiving the position and the magnitude of the 
current element at another point of time when it is judged in the step SP5 that the end 
of the excitation current display is not instructed (step SP 6). Then, the step SP6 goes 
back to the step SP2. When it is, on the other hand, judged in the step SP5 that the 
end of the excitation current display is not instructed, the process goes to END. 
[0073] 

As will be seen from the foregoing description, it is to be understood that the 
present embodiment of the organism magnetic field measured outcome display method 
according to the present invention can display the sphere having a diameter 
corresponding to the magnitude of each of a plurality of the current elements along 
with the semi-transparent organism structure reference image for each of 
corresponding points on the time series, thereby making it possible for an operator to 
easily recognize the three-dimensional propagation direction of the excitation current. 
[0074] 

[Tenth Embodiment] 

Referring to the drawings, in particular to FIG. 23 of a block diagram, there is 
shown a tenth preferred embodiment of the organism magnetic field measured outcome 
displaying apparatus according to this present invention. The organism magnetic 
field measured outcome displaying apparatus is shown in FIG. 23 as comprising: an 
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analyzed outcome storing unit 600 for storing therein current elements analyzed 
outcomes, a circle assigning unit 700 for assigning a circle having a diameter 
corresponding to the magnitude of each of current elements to each of the current 
elements, a tomogram image displaying unit 800 for displaying an organism tomogram 
image, a circle displaying unit 900 for displaying the assigned circle having a diameter 
corresponding to the magnitude of each of a plurality of the current elements along 
with the organism tomogram image. 
[0075] 

As will be seen from the foregoing description, it is to be understood that the 
present embodiment of the organism magnetic field measured outcome displaying 
apparatus thus constructed can display the circle having a diameter corresponding to 
the magnitude of each of a plurality of the current elements along with the organism 
tomogram image, thereby remarkably increase the visibility of the excitation currents. 
Furthermore, the present embodiment of the present embodiment of the organism 
magnetic field measured outcome displaying apparatus can display the circle along 
with the organism tomogram image on the time series, thereby making it possible for 
an operator to easily recognize the three-dimensional propagation direction of the 
excitation current. In the previous seventh to tenth embodiment according to the 
present invention, the circle may be replaced with a square, a cube, or the like as long 
as degrease in visibility is admitted to some extend. It will be understood by those 
skilled in the art that the foregoing description is in terms of preferred embodiment 
of the present invention wherein various changes and modification may be made 
without departing from the spirit and scope of the invention, as set forth in the 
appended claims. 
[0076] 

[Effects of the Invention] 

As will be seen from the foregoing description, it is to be understood, the 
tomogram image processing method as set forth in claim 1 comprising the steps of: a) 
carrying out an interpolation operation on a plurality of tomogram images of a test boy 
sliced in an area at predetermined intervals to obtain stereoscopic image data elements; 
b) dividing said area into small stereoscopic segments; c) counting the number of said 
stereoscopic image data elements for each of said small stereoscopic segments; and d) 
binarizing said small stereoscopic segments on the basis of the number of said 
stereoscopic image data elements, has an effect of obtaining the stereoscopic image 
data elements corresponding to all of the small stereoscopic segments by carrying out 
the interpolation operation on a plurality of tomogram images of the test body even 
through some tomogram images sliced at the intervals should fail to be provided. 
Further, the tomogram image processing method thus constructed can obtain a 
simplified stereoscopic image by the reason that the small stereoscopic segments are 
binarized on the basis of the number of the stereoscopic image data elements. This 
leads to the fact that the stereoscopic image obtained and displayed by the tomogram 
image processing method as set forth in claim 1 makes it possible for an operator to 
easily distinguish an area and/or a stereoscopic segment where an internal organ or the 
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like is placed. 
[0077] 

A tomogram image processing apparatus as set forth in claim 2, comprising: 
interpolation carrying out means (3) for carrying out an interpolation operation on a 
plurality of tomogram images of a test boy sliced in an area at predetermined intervals 
to obtain stereoscopic image data elements; area dividing means (7) for dividing said 
area into small stereoscopic segments; stereoscopic image counting means (8) for 
counting the number of said stereoscopic image data elements for each of said small 
stereoscopic segments; and binarizing means (9, 10) for binarizing said small 
stereoscopic segments on the basis of the number of said stereoscopic image data 
elements has an effect of obtaining the stereoscopic image data elements corresponding 
to all of the small stereoscopic segments by carrying out the interpolation operation on 
a plurality of tomogram images of the test body even through some tomogram images 
sliced at the intervals should fail to be provided. Further, the tomogram image 
processing apparatus thus constructed can obtain a simplified stereoscopic image by 
the reason that the small stereoscopic segments are binarized on the basis of the 
number of the stereoscopic image data elements. This leads to the fact that the 
stereoscopic image obtained and displayed by the tomogram image processing 
apparatus as set forth in claim 2 makes it possible for an operator to easily distinguish 
an area and/or a stereoscopic segment where an internal organ or the like is placed. 
[0078] 

The organism magnetic field measuring method as set forth in claim 3, has an 
effect of obtaining the data on the magnetic field source for each of a plurality of 
vertexes of said small stereoscopic segments by the reason that data on said magnetic 
field sources included in the results of the Biot-Savart law calculations is corrected 
until the difference between said calculated magnetic field value and a measured 
magnetic field value becomes vanishingly small. Furthermore, the number of 
unknown factors to be estimated can be largely reduced in comparison with the 
conventional method by the reason that the position of each of the vertexes of the small 
stereoscopic segments is known. The organism magnetic field measuring method 
thus constructed degreases the number of measuring points, thereby as a whole 
reducing the operation load. The organism magnetic field measuring method as set 
forth in claim 4 has an effect of obtaining accurate information on the magnetic field 
source by the reason that even magnetic field source vectors assigned in positions 
adjoining to each other, equal to each other in magnitude but opposing to each other in 
direction can be calibrated. Furthermore, the organism magnetic field measuring 
method thus constructed degreases the number of measuring points by the reason that 
one vertex of one small stereoscopic segment serves as the vertex of the adjacent small 
stereoscopic segment, thereby as a whole reducing the operation load. 
[0079] 

The organism magnetic field measuring apparatus as set forth in claim 5, has 
an effect of obtaining the data on the magnetic field source for each of a plurality of 
vertexes of said small stereoscopic segments by the reason that data on said magnetic 



30 



field sources included in the results of the Biot-Savart law calculations is corrected 
until the difference between said calculated magnetic field value and a measured 
magnetic field value becomes vanishingly small. Furthermore, the number of 
unknown factors to be estimated can be largely reduced in comparison with the 
conventional apparatus by the reason that the position of each of the vertexes of the 
small stereoscopic segments is known. The organism magnetic field measuring 
apparatus thus constructed degreases the number of measuring points, thereby as a 
whole reducing the operation load. The organism magnetic field measuring apparatus 
as set forth in claim 6 has an effect of obtaining accurate information on the magnetic 
field source by the reason that even magnetic field source vectors assigned in positions 
adjoining to each other, equal to each other in magnitude but opposing to each other in 
direction can be calibrated. Furthermore, the organism magnetic field measuring 
apparatus thus constructed degreases the number of measuring points by the reason 
that one vertex of one small stereoscopic segment serves as the vertex of the adjacent 
small stereoscopic segment, thereby as a whole reducing the operation load. 
[0080] 

The organism magnetic field-measured outcome display method as set forth in 
claim 7 of displaying an organism magnetic field-measured outcome along with an 
organism structure reference image, has an effect of making it possible for an operator 
to recognize, for example, the distributed condition of the maximal organism magnetic 
field-measured outcome in comparison with the case that the magnetic field measured 
outcome is displayed with the direction and the length of the arrows. Moreover, the 
organism magnetic field-measured outcome display method thus constructed can 
dynamically display the organism magnetic field measured outcome, thereby making it 
possible for an operator to easily recognize, for example, the time-variant fluctuation 
of the distributed condition of the maximal organism magnetic field measured outcome. 
The organism magnetic field-measured outcome display method as set forth in claim 8, 
in which the reference graphic is in the form of a circular shape, has an effect of 
considerably enlarging a clearance between the reference graphics, thereby making it 
possible for an operator to recognize the organism structure reference image as a 
background. 
[0081] 

The organism magnetic field-measured outcome display apparatus as set forth 
in claim 9 for displaying an organism magnetic field-measured outcome along with an 
organism structure reference image, has an effect of making it possible for an operator 
to recognize, for example, the distributed condition of the maximal organism magnetic 
field-measured outcome in comparison with the case that the magnetic field measured 
outcome is displayed with the direction and the length of the arrows. Moreover, the 
organism magnetic field-measured outcome display apparatus thus constructed can 
dynamically display the organism magnetic field measured outcome, thereby making it 
possible for an operator to easily recognize, for example, the time-variant fluctuation 
of the distributed condition of the maximal organism magnetic field measured 
outcome. 
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[0082] 

The organism magnetic field-measured outcome display apparatus as set forth 
in claim 10, in which the reference graphic is in the form of a circular shape, has an 
effect of considerably enlarging a clearance between the reference graphics, thereby 
making it possible for an operator to recognize the organism structure reference image 
as a background. 

[Brief Description of the Drawings] 

The objects, features, and advantages of an imaging device correcting 
apparatus according to the present invention will become apparent as the description 
proceeds when taken in conjunction with the accompanying drawings. 
[0083] 

FIG. 1 is a flowchart showing a preferred embodiment of the tomogram image 
processing method according the present invention. 
[0084] 

FIG. 2 is a schematic diagram showing an example of a tomogram image. 

[0085] 

FIG. 3 is a plan view diagrammatically explaining a relationship between 
small stereoscopic segments belonging to a small stereoscopic segment group to be 
displayed and a small stereoscopic segment group to be displayed. 
[0086] 
[FIG. 4] 

It is a perspective view corresponding to FIG. 3. 

[0087] 

FIG. 5 is a block diagram showing a preferred embodiment of the tomogram 
image processing apparatus according to the present invention. 
[0088] 

FIG. 6 is a block diagram showing a preferred embodiment of the organism 
magnetic field measuring apparatus according to the present invention. 
[0089] 

FIG. 7 is a block diagram showing another preferred embodiment the 
organism magnetic field measuring apparatus according to the present invention. 
[0090] 

FIG. 8 is a schematic diagram showing the detailed construction of the part 
corresponding to one hierarchical type perceptron. 
[0091] 

FIG. 9 is a block diagram showing another preferred embodiment the 
organism magnetic field measuring method according to the present invention. 
[0092] 

FIG. 10 is a block diagram showing yet another embodiment of the magnetic 
field source measuring apparatus according to the present invention. 
[0093] 

FIG. 11 is a schematic diagram explaining a process for extracting an adjacent 
vertex, for which an estimated vector component of a current element is to be 
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distributed. 
[0094] 

FIG. 12 is a schematic diagram explaining a process for estimating a current 
element carried out during a redistribution operation. 
[0095] 

FIG. 13 is an illustration representing information on vertexes with depth of 6 
cm obtained based on MRI pictorial images. 
[0096] 

FIG. 14 is an illustration showing an estimating outcome obtained at the time 
35 msec prior to the occurrence of R wave. 
[0097] 

FIG. 15 is an illustration showing an estimating outcome obtained at the time 
33 msec prior to the occurrence of R wave. 
[0098] 

FIG. 16 is an illustration showing an estimating outcome obtained at the time 
30 msec prior to the occurrence of R wave. 
[0099] 

FIG. 17 is an illustration showing an estimating outcome obtained at the time 
25 msec prior to the occurrence of R wave. 
[0100] 

FIG. 18 is an illustration showing an estimating outcome obtained at the time 
20 msec prior to the occurrence of R wave. 
[0101] 

FIG. 19 is a flowchart showing a preferred embodiment of the organism 
magnetic field-measured outcome display method according to the present invention. 
[0102] 

FIG. 20 is a flowchart showing another preferred embodiment of the organism 
magnetic field-measured outcome display method according to the present invention. 
[0103] 

FIG. 21 is an illustration showing examples of excitation currents displayed 
by the embodiment shown in FIG. 20. 
[0104] 

FIG. 22 is a flowchart showing another preferred embodiment of the organism 
magnetic field-measured outcome display method according to the present invention. 
[0105] 

FIG. 23 is a block diagram showing an excitation current displaying apparatus 
as a preferred embodiment of the organism magnetic field-measured outcome 
displaying unit according to the present invention. 
[0106] 

[Description of Notations] 

3 Tomogram image interpolating unit 

7 stereoscopic segment assigning unit 

8 pixel number counting unit 
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9 judging unit 

10 display-stereoscopic segment assigning unit 
20,200 adding unit 

30,300 error calculator 

41i correction quantity calculating unit 

42i vertex extracting unit 

43i standardizing unit 

44i redistributing ratio calculating unit 

45i redistributing unit 

46i estimated value updating unit 

40, 500 data collecting unit 

101,102, ... ,1 0m hierarchical type perceptions 

111,112,... ,11m Biot-S avart arithmetical calculating unit 

llla,112a, ... ,llma, 121,122, ... ,12m correcting unit 

700 circle assigning unit 

800 tomogram image displaying unit 

900 circle displaying unit 

ABSTRACT OF THE DISCLOSURE 

PURPOSE: 

To obtain a solid model from tomogram images 
CONSTITUTION : 

The construction of the foregoing methods according to the present invention 
comprising the steps of: performing interpolation arithmetical calculation on tomogram 
images taken at respective predetermined intervals of an object such as, for example, 
internal organs to obtain stereoscopic image data elements, generating small 
stereoscopic segments in correspondence with the tomogram images and/or 
interpolation tomogram images to count pixel numbers contained in each of the small 
stereoscopic segments, assigning the small stereoscopic segments to a 
display-stereoscopic segment only when the counted value of the pixel numbers is 
larger in number than a predetermined threshold, and visually displaying only the small 
stereoscopic segment assigned to the display-stereoscopic segment. 
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METHOD, LIVING BODY MAGNETIC HELD MEASURING RESULT DISPLAY METHOD AND DEVICES 
THEREFOR 

(57)Abstract: 

PURPOSE: To obtain delicate image data, based on a tomographic 
image at every prescribed interval by dividing a space which 
becomes a processing object into prescribed small stereo scopic 
areas, counting stereoscopic image data contained in each small 
stereoscopic area, and binarizing the small stereoscopic area based 
on a counting value. 

CONSTITUTION: Based on a tomographic image at a prescribed 
interval, an interpolating operation is executed and stereoscopic 
image data is obtained (SP4). A space which becomes a processing 
object is divided into prescribed small stereoscopic areas, 
stereoscopic image data contained in each stereoscopic area is 
counted (SP5), and based on a counting value, the small 
stereoscopic area is binarized. Even if only the tomographic image 
at every prescribed interval is obtained, the stereoscopic image data 
is obtained so as to correspond to all spaces by executing an 
interpolating operation, and also, by binarizing the small 
stereoscopic area based on the number of stereoscopic image data 
contained in the divided small stereoscopic areas, a simplified 
stereoscopic image can be obtained. Accordingly, in the case where 
the displayed is executed visibly, the part corresponding to internal organs, etc., and the part corresponding 
to the space can be recognized simply. 
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*ffftoTiCft«r-**^S»f«flMfH#® (3) 

£hbk##& (7) K^«nfe*/jNi{«ii«t^s 
(8) tHRatJStJv>T/hitft9i«*2ffl[{t-rs2 

fifb¥© (9) (l o) ^^#t?ci:^!if®i:1-§»fii 

Mia 

[Mr 3 ] m-m. i ® K*«fflSL&i£fc ± t> # 5 n 

[«t*S4] *fflbfeSfc:S^<fflliEfcftAT, #/h 

[w*s 5 ] mm i o«f*««ra^i£t * 5 # 6 n 

a^K^S (111) (112) • • • (11m) (1 
0 1) (1 0 2) • • • (1 0m) ts &m.W^Wi (1 
11) (112) • • • (11m) (10 1) (10 

2) • • • (iom) ^6ta*sn*adH8**jRatta 

©ir;si»}¥#® (20) (200) mmtm^ 

& (2 0) (2 0 0) fr5ffl£SftSJR»faJBS»fc« 

(3 0) (300) t, WffilRMfc»-3v"»T&iJW#® 

(111) (112) • • • (11m) (10 1) (1 
0 2) • • • (10m) lCi5if&m&miim<DffiE%fi 
&5*iIE#I£ (Ilia) (112 a) • • • (11m 
a) (12 1) (12 2) • • • (1 2m) t, fifiE# 
IS (1 1 1 a) (1 1 2 a) • • • (1 lma) (1 2 
1) (12 2) • • • (12m) fciSfluEtfSSSttfc 

(4 o) (5 0 0) h^tsct^mtt 



o#jsai«4)-o«iE«**ffl-rsflEjE*»m#a (4 1 

m&mm%BiEtt»m£Mtii^m (42 n «• 

(4 3 1 ) (4 4 i) (4 51) (4 6 i) t^ttffe 
[M*«7] 4#aSij|II£ilS««4ft«Bfi#j!S»i:« 

® x*sas u ^mmmimm^ % 1 1 & 
20 ©^-r x*w&r s *&h«im xsmmm (700) 

^«3i#M«^S^-r§^ft«lS^#@ (8 0 

0) t, ^.wm^mimm-mmim^nrc^x 
o*ffiHjg%as?i-s*ffiHjga^#® (900) fcs 

[ft*J 1 0 ] *ffiH^R-ea&«»*S 9 £IS«cD 
[fS^lfffl&iKBJ!] 

[0 0 0 1] 

a^-r § * <r>m M%±mmmfe%£$km*j5m & * i? 
[0002] 

«!8?«St«ft»lc, MR I Bff ^©WMff ^RlSKl 

«s 3 ^ozmzm.^ l, 5i«jis^ 1 3 Xftrnm^t % 

[0003] tfc, 4ft^iJ^aiJS1-S^Si:L 

50 ^\znm^r^m\^m^^xsmmm.-hmwmMm 
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fcj:iraiRitt5*ifia-rs#i£-efc5o 

m\t. mmmmtfr (asm 3 o~i o o@g©«^t 

[0004] tfc, s#««fKiiis*ts-5<aa»sjB 10 
■^<3>*-Ty^a^a<»i5nTv«So ^lt, is* 

jB*Tfc/flV>5*£*:»* MR I HfiLhfcMflE'** Wl/fc 
[0 0 0 5] 

[^Bj^S&LcfcSt-tSHg] ±EJ8SaS!a^SO 

(WfciMT*^iHioffl«ffiii) ff£<fflB?#&vMzre 
(Dtitm*mmti:imt& c tmmic ntsotu 

[0 0 0 6] 3 ^WcS^" 5 73 Soffit ft 

^KWffl-rs^i&oai^ffij&sffiL^Ji^Ktt, as 40 

[0007] tit. mmmmmmffi.mfr?%z>£iixM 
Wt, 3 0-1 0 0 ggowiautsss-r s c # & 



air-rtiffx 2 0 0~1 0 0 O^+iM/gSOTW-fe 

[0008] *fc, ±i3Tn-v-y7 p a^ Hiing* 

HPS, SL^itcDMRJt^^ Wl/g 
MR IHfltetf^-^yT^atflfcfcBlfcjiTfc, 

[0009] 3>^-vy^a^*»«a^fta»si3£ 

[0 0 10] 

mm<o mi c o^tt±iB©iHiiijat**Ta«n 
[0011] 

[0012] M^2<DKrAtt^&s^ mfemmm 



(4) 



#M¥ 6-28455 



[0 0 13] ±IBil2©gW£«jjjc-t5fcJ6©, !t*« 
3 0£f*?^»|j££i£tt» flf&g 1 OKMfiWBBl^afc 
J: D » &ftfc/Jvfc{ra#©Kja®ffim* skis^tii 

*Jt3¥ffi Lfcj&ifff ifJfE i: ®g£ ifffi U HtB Lit 
[0 0 14] »3^4©£{ttK*ai3£75f}£H\ SffiLfu 

[0015] imm5<D&mmm^.mmii. hj&su 
(Dmmmmwxmc «t 0 6 ftfe/hftfl^aoiaaoffi 
miisiES-jv^TaasafiiSfflfflitstJ < as® 

*frft 5 «uE#® £ , ffilE^g fc «t 3 fiuEtfSSS ttfc*S 
[0 0 16] W*S6<D£«raSi!EW^ fluE#® 

at, #j«g^©«iE«ta^v^T, nznmmft 

^afcfc^trfcOfcfflWtVSo 
[0 0 17] ±fE83®I«£38jS-r*fc«><s^ 

[0 0 18] «5f<]19 04f«^iJ^m^gBa, 



[0 0 19] 

»tscyv^T/jNitft«i«« 2 fiifc-r * c fc t * t> «fc 

£ftfc:fc#*fc8Si: LfttfoT, nffflW 

[0 0 2 0] ft, ±fBSW»Wi:LTB, tttttflSft 

[0021] «*3S2o(ff*«sp^ai?a&n«s m% 

LT*3<o fit, Aft«T-*tHR#atJ:t)» E# 
so £ ft fc&/J vfcttfflWfc # £ ft 3 £ x- * *fHR U 

2 mit^mz j: d » ttsffiica^v^T/j^ft^^ 2 fi 
o3fe(ta»si£7jffi -efeft^ i o»f»« mmz 

ffifc i 0 # 5 ft^fift^oiljSioffi«fll«fc«-3v^ 

uoaw^ff a v \ «»dns«« jr%!iq3* t t# 5 ft § 
ufcMfcas-^^T^awsic^sftswssfflfflisffiiE 
T^6#«fl^t^sft«i»w«fflfiB*£f*aa»iiisis 

ffi^^ft § S T?»iE*ff ft a C £ «t D s #/J^fSS 

[0022] »^S4 0^ft«Jiais^i£-efefttf, § 



(5) 



ft HB¥ 6-2 8 4 5 5 



co o 2 3] »*S5©^#i8»aissn?*n^ ft 

or, sswi^acj;^ ^idw^gj^sm^sns 

**ns*»bi^«taHifHIIffl[i:K:as^->TBl!IJ|[ 

5 o SiE^at «fc SffiIEiQJi£&gEI»7£«-?T&o fc« 
K, ffiiEISftiDafc^atJ:?), ffiGE^afciSffijEsWM 

co o 2 4] »«S6o«»ssii®einffen«, *hie 
n a *> c 1 1 j; 0 , *m l fcs t so* < siESff a a © 

■T SKjSil1B-PSE»*fT 4 5 c £ lz «fc t> ffi££ff & 3 o 
CO 0 2 5] £ftS&if$iJ£#i££o^T2££i¥fflfc!M 

co 0 2 6] 1 o©/Jnsm«o 1 -ocDmmz, 



CO 0 2 7] S5fc$fc, «»«TftV>]IjSfcfetJ5« 

§0 3 ^ 6 <omut±m<mfiK.&^ 

Tft?ftfcfc©T?& !), SI£l»&Hf j£fc:feV->TH:*fc]0 

f « <: i: «fc ?> 5^o«3(EJRM-<0WSIfioJtS%3i^-p 
JtS««**«:jK*iai£i|g*kLTSWTf#5„ 
20 C 0 0 2 8 ] flt$9 7 ©^^^SiJ^Sg^STfe 

^rs*&tjt«LT, mar, «*©4ftawiaisjg 
co 0 2 9] »*B8©4ftiK«ai3eig»a^js-efe 

5, *ffiH^lRl±OHilH^S:»)**<T?t, ffiRfcft 

§ i: t «> ^fia^a^^afc ± d ^mmmm 



[0 0 3 0] ttSa 1 0O£#a»jlI£!l8«S^Sfli? 

fthtf, mass 9 ©Austin*., ^mm^^h^ 
s £#*jS#jr* ©teg jwsjb k a 5 o 

[0 0 3 1] 

1>T, 5 nfcW««fc*3 < SP^S^T So Tffi 
HSfJlff^f, Xr»y7SP3fc}3V>T, #A5tifc»f 

fc*FfcU *xy7SP4fc*sv>T, *?{fc£ftfcIE 

7^©^l^*/d£ U fro 1 atflE&JB® 1 Ht^L 

ffl»f»««*^1-5IiaiSO»SIHRL, 7f7^SP6 

WBXfy^S P7lcfcVTl^Jni.'<t]iSft 

7S P 6 fcfe^TtHSffitfHffii 5 t>±£ < ftv^Wgij 
*nfe3t&ftfcov>TfiXr«yys P 8t&^TSS£ 

nftv^ft cwt* i^^si^fcftrrs) tilt) a 

ftt>tifeftfi» Xf-yys P 9 fcfcv>T£T©:4£f*fc 
'^SPl Ofcfc^TffiOjHfttfcaftU nvxrv 

^sP5offla*frfta. tfc ^ysP9t^ 

[0 0 3 2] fit, 4##fc£S:ft5IIij|5»fcft*.T, 

iw (D±mmm t *rr § »fSfiH«so»J^te ± 
H 2 e. H 4 t & £ KSftsaa^iS * $ £ 1 K 
o^L^^nfcgP^A^OUgs^LTV^o H2 



(6) #M¥6-2 8 4 5 5 
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* V ^fr&fc {± H 3 1 Tffi*f L T^f <fc 5 tcSzttftW £ 

*a^irtriti:K«i?)ST5nSo h 4 3 

[0 0 3 3] ffi, CO-aiO«BlS±TO»fM«feJ:tf 
>y F-^xVl^LTW^cg^tSo flu flttlS 

[00 3 4] 

[H»j2] m5itc(Dmm<Dmmimmmmco-'mm 

mm&m 1 Bfjift{g&9 1 ^ sfflSflafcasK 

m\mmiirz>m 1 KS{faK^2 1 % aK^nfcifs 

< ffiEWaafcffft o TffiHSfff 

t , mmm mm® z ^xfKm«mi» 3 n&gm&ie 
\mm.mfti£&&m 1 Rwmm 5 ^ 1 stsj»gp 
30 5\z^mm\mm^miLT^h\c^Lx, mm 
mmm 1 $^(iMH0r«ft^i!4 t^mmmm * 

*SLT3£&ft*ffl5aT*i^(ttfl0 3Ti5 7i:, sou 
«5 1 0 fc, ■3R»^m3£OBBfil«fe u < 

swsitssfiifipas 12^ n 2 sffsft a^gp 6 , 

H«ti-agP8, wsygp9, s^ii^sijagp 

»ff«-&SS2Sa#J»i5l 3t, ^2g«SlJ«l 3 
50 [ 0 0 3 5 ] Lltff? T, 0 5 08j£0BrWS&3&B 



(7) 



<8WF6-2 8 4 5 5 



[0 0 3 6] 

[HSl« 3 ] H 6 tt c ©fg0^©4fr^^ffl^g*O-^ 

im%f&t7uy9mTi&*). mmm<D, £* 

;l/®i£B0©»W*fTa5 If* • <y h l 

11, 112, • • •, 1 1 rat. if* • 
a-yMii, 112, • • •, 1 lmfretHft^ft 

>jh? - a- >y h 2 Ofr5fflft^5S«flnJHg*0 j 

(t) tmm'W->tLT<DMmmms j co t 

*Xh t LTPi#©S*»tH-r SMSSORS 3 0 fc, w 

mzftrcmm^T, tr* • -y-^-^ga- y h i 

11, 112, • • •, 1 1 mfcfeVTJtS^nTVS 

maauT-o^^h^^ffijE-rssjEapi 1 1 a, i 

12a, • • • , llmat, W - ■y-A-^grL- 
4111, 112, • • •, 1 lrat^TiSSn 

Ttttal-SflNBiRifci-'y h 4 0 tZHLT^Zo 20 
±IBlf* • Mil, 112, • • 

•, 1 lmli, RfSflu ^10BfS«M^SfeL<ti 

0 5 vmmwmsmtc «t t> » 6 nfc» n^idf mew 

BteStifcC fcfcj£«FLT, K»M*fc»3v>T&lf* 
• MdS££ttTV>3, If* • If ^ 

-;i/offiai<oadi*ff ft 5 1 1 1 gtsgggg 30^ 
sm^snsfiuaRMd j (t) {=sj co -oj 
(t) } m^tirccticm^Lzm.w^^n^ 
mmfrv^ h frtitftztmfemmtf'p* < ft § ± 5 \z 30 

SlE-TSo £7^ If* • ^-MOH^y Mil, 

ii2, • • •, i lmttwenrawwfciWFrsj: 

[0 0 3 7] ±f2©§$©£«^$J£SB©ftffl&# 
OtfeDTfeSo SWTfcffcfcftSfiBfOJ (t) «\ B# 
SO t fc&/h££#<Bl&£W> 3 ^Ttlffigffi x i , y i , z 
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1 (itiiEosas, othu) ts&xsmmfrwv 

M^Pxi, Pyi, PziWJigi (If*- 

[0 0 3 8] bfctf ot\ BfgiJ t »J:tf#/jN3t^{tOS- 
]ijt0 3^JS««£mffl<73lf* • f/HMifflti-y 
Mil, 112, • • • , 11 mfcSa&LT^ft^tt 
ilftgl. g2. • • •. gmOaWSfrft-sTHSfi 

2 o fcflaei-sc i:K i t) JRSEtaJMO j ( t ) 
C£tfT-t3o fit, aWfi^maSRfrO^^ 

JWIO J ( t ) (ilgR&DtHgMI S j ( t ) t ttSft So 
LfctfoT* KS»JH83 OtfeVT^SSEOfHIIMiiS j 

(t) tjnainiiwBo j (t) tvmwmu wm 
tiK&zmfemm d j c o t i/nf* • 

any Mil, 112, • • • , 11 mOffilEgP 1 1 
la, 112a, • • •, 1 lma(C7^-K/W^ 
U d j ( t ) *VJn$ < ft£ i a fc£lf * • +)" 

[oo3 9] ±E-iiosaa*sa-rntf«iSRMd j 

( t ) # < ft o , ov>t imjmm d j c t ) #s 

(£0£ft5©T\ COH&^fcfcVTlf* • im-;]/** 

«t 0 HX* bTtitfrr * c i: t i. D £/h3Kfft®#M;: 

ft. S£ilMW«i&Ej (t) *&£T?£»*-nfcr, 

E j (t) = (1/2) {S j (t) -Oj (t) } * 
[0 0 4 0] 
B&l] 

[ S j ( t ) - 0 j ( t ) ) 



[0 0 4 1] ZLT. £lf* • tf/WMSfflti-y hfc fit, e ktt^^ h;l/fig^a i kO^S^V (ffiiEy , 
fc'tt5^4;l/^OffiIE^«f^T?£t»^^Tffft 40 T*h%o 

*>ct t?m, , m^mmmmmm^m'hicr^^ [0042] 
hmftvm&m2 tm^zfitts c t^ts 0 as 2] 

a i k-a i k-Ek 19 E j (t) /9 a i kj 
-aik-ek (OEj Ct)/30j (t)] (30j ( t ) / 
Jalk] - aik + ek (Sj (t)-Oj ( t ) ) {90j Ct) 
/9a Ul 



[00 4 3] f fc, liiDllO j (t)Oaik£J: [0 0 4 4] 

sffltt5MBa» 3 ? # 6 n§ ot, a 2 as 4 mm imi 

T*f§o 50 
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30 j (t)/3aik = 3 (Eg i (t, 
• "aiL)| /9aik 
- d Ib i (t, ail. a L 2, a i 3. 
9a i k 



ail, ai2, a i 3, 
• • • a i L) } / 



iik- 



[»4] 

ik + Efc |Sj (t) -0j (t)} id (g i (t, 



i 2, 



i i L) j /a 



1 k] 



[0 0 4 6] LfctfoT* ft4©jZra*frft5Cfct«t 

fiSc^#§ci:^-et§o ^ JtSKMd j (t) ©fi 
t?tfIE0»&fcH:«uEffi*£fc, SSKMdj (t) © 

[0 0 4 7] $fc, H6 0SBfc2S-5<J!BlttHl*fc 
tt & D Sv«*?ttO<IxK*M)IMi%*tiW5 O 

[0048] set, a»ittaifi«»sfe«>oi!8£jsas 

[0 0 4 9] il©^ac^V^T?F,tp»fcSJi§-rSc, $ 40 
5> SIISjSi£»i*&ffSavvjNi^ftfcfe»S««E3R so 



W-ao^^h^fciRSL/TUSactfcftSo fit, 
0^^h;l/tiRjRLfe«SKJ(yttoV'>Tti««Oifl®a 
SiS^S t ftS COTS 6 £>, M«'>SV^gjiCCf 

[0 0 5 0] 

^SSM^-T^n-y^a B8ttlocD«»S!/*— fe7 
^X* l tH*©»»fflOPg*S^-- tyhnvi o 

p ( P =i, 2. • • •, m) t, mmm^-tfvu 
y i o pfrs®BBSffltii;&g i j snainiw-rss/^ 

• a->y h 2 0 0 t, S/^V • ary h 2 0 0^6 tH* 

*ns*aiinjHg*o j (t) mm^z-ytLTCD 
mms j (t) t*xtitLT-ffi%<Dmmiiitzm 
9kmmmmnm3 oot. ^mmm^~^vuy\ o 

5S»#«lfcBg|$ 1 Opr (r = l, 2, • • • n) 
SKA®/*-- \i7Yuy\ 0pOffl«iHKfffli5i 

0 p r^SOffi^^Sffl^^fc^D^ffffl^nT 

M#OS*Jlffl-r«fiMK^MIRS«»«3 0 p r t, M 

mmmmm3 o ofc£ DWffi£n3M££t«# 
mm&m3 o P r ic^^mm^^mizm^rm 
m^mmm^-^fvuy 1 o pfcarata-Hso^s 
*ffftfe*«ffiia»H»^siJ6 opt, mmmmm 

WS3 0 0fcJ:t)J(tiJSnfeliSHinJWSSO j (t) t 
LTOtHfflfflSj (t) i:0gfe,fctfffi 
fifcM»tHW l 0 p r fc«fc !>*ffiSnfcffl«$Mtfc»-3 
^TBJiS!^— fe^ hP> 1 0 pOAWifcfeW-SArt 

fcSiiE-rsttKEffls 1 2 p t, mm^ 1 2 P tiamttmm 
^§§13 6 0 P fc*aft-r* a 1 1 k % aftsn^ffiiEiJ 

1 2 p*fcB«tt#Hft£S»6 0 p^3fsma (IS* 

-t-^t/jN* < t? t s nso ^ttsa»f¥?-&s«o»» 4 

0 0t, SlEgP 1 2 p fc«t jE^BfitiaR^ga^ 

ftfcc tzmt Lxmmm^-^fhuy 1 0 P © 

LTffi^-rS'ltSiKfta- y h 5 0 0 t^LT^ 

§0 ±m&mmm' <-tyhn^iopii, to 



(9) 
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Coo5i] sfc, anK3-H8&3S6 o P a, m 
ar, mmmmmmwm 3 o o t j: & juasnssfc*-? 

<Pf»S^— fe7hn>l OpO^S 

yn;^y— > 3 >^S) &J:tfflraKHfB@Kffgg3 0 

p r t £ D WlUSftSgK:S-3< B*3!/*— fe7 h n v 

[0052] set, «ie»i 2 pti^n^nascffi 
LTfe t) „ mmmmmm 3 o o k ± t> »tn 
snfcMfeitfffla^Horttjas 1 0 P kfcionasn 
feM%x*fcbTKs-r«5fesia»*aiE-r*o ±e© 

PgHH^— b7°hny 1 o ptfevre* • •9vv-;KD 
McSo* < ^S#s57 L TV S <Z> ? & S 6 , *£n<D 

ffiWfflap 1 0 p r fc«fc D»*S1"Sffla5MI*»S 0 i& 

c ottUt fc^Tfel W 4 0 fc J: D OMIKMHR^Sas 6 
0 pfcgftU ^LfcOI«#HS3*B15 6 o p^ms 

H»fe *tf c ©S8»<Dffl«^HaoifW2:«ilia«S* 30 
3i^7*5J:5h:#Pi»l!^— b^hnyi 0p« 

[0053] ^ nm^mmm^— t-fhuyi o P 

£*«fcJE«T**So fiLhO«fc 5 1 LT&ga#e#S 
7Lfc» fitX#HS&B$6 0 PfcfW.Tft9flS4 
0 0 1 «t t) Jl IE8& 1 2 p %mR* S k fe KflBBJR*a 40 

-■^500 WtS^nH: £ < % m 6 0^jt«9i:RI« 

[0 0 5 4] 

im<Dmmm^tytt 9 ^h?s t> > a 7 & a 8 

— b^hnyft£«t3-6nTVSfi«5HtJHB»i 0 p 
r IZttTLT, 1 Art®#£ffeMtf£ttg{bS-&S A^g 

im 1 5 p t, i A*os{tw«fcfeftsffl**ffiRp-r 
sm^fflfiris 1 6 p t, nmiKD^w-mzm-ftm so 
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i 7 p£Kttfcja<z>*"Pfc5o 

[0055] Lfttfot, coHfifflBio^fca, 

[0 0 5 6] 

custff J 5 ] a i o it c <D5m<DmMw$mm<D 2 6 

#*^bTVS„ if* • ■9-^-;l/OffiaiOifl[ 

• Vrt-AM&fr-y Mil fcfcVTtiSSftTVS 
WORK-®"** WVj£#*^LT«^£fcLTfflrt 

■T5«a0IUT-I|5»3.^y h40a^ HI <Z>£}£i:fctt 

0 5 OSHK <fc t) » ?»nfc 3 ^TtMItt^S^ LTif 
iEMAx, Ay, A z£jrffl"rSflfiE**mff4 1 i 

»ffl5nfc«iE»Ax, Ay, AzOKWKiaS^ 

^fc^sssn^j4ffi-r§]i^affia54 2 1 1 % ^je 

IAx, Ay, AzK3S"3VT£SSiaEiE^^h;l/*iE 

mct^wm^ 3 1 Mb^nfcfiiE^^ 
^it*wma5 4 4 1 jeesnTv^maasfro^^ 

»KJi^ffiO|IE»SffftdSE^4 5 i BE^ 

354 5 i tc£Df#E#2nfcBE^W*^VTBSir 

ifg|54 6 i i:, ^TOMc-OV^TffiEMA x, Ay, 
A z AWt/jN? < £5 fj-pe* • 

mm, immw-mA 1 1 , n^affi^4 21,1 
mm a 3 1 % sE#jfc**tH»4 4 i , se^4 5 

1 fe<to : «sfiMffgi54 6 i osas^sa^^ssasij 
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